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ABSTRACT 
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Memorandum.  The  conduct  of  the  DRT  (used  for  quantifying  speech 
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INTRODUCTION 


1.1  The  Diagnostic  Rhyme  Test 

The  Diagnostic  Rhyme  Test  (DRT)  is  used  extensively  for  assessing  the 
intelligibility  of  military  communications  systems  and  has  become  an 
accepted  NATO  standard  for  testing  Linear  Predictive  Coders.  A 
discussion  of  the  acoustic  and  phonetic  research  on  which  the  test  is 
based  is  beyond  the  scope  of  this  report,  conssquently  only  a  brief 
outline  of  the  test  is  given  here. 

The  DRT  voceQjulary  comprises  ninety-six  minimally  contrasting  rhyming 
word  pairs,  the  initial  consonants  of  which  differ  only  by  a  single 
acoustic  feature,  or  attribute.  There  are  six  such  attributes; 
Voicing,  Nasality,  Sustention,  Sibilation,  Graveness  and  Compactness. 
As  an  example,  the  attribute  voicing  is  present  when  the  vocal  cords 
are  excited;  in  the  word  pair  "veal-feel",  the  consonant  "v"  is  voiced, 
but  the  consonant  "f"  is  unvoiced. 

The  DRT  is  implemented  using  the  following  procedure.  The  192  word 
vocabulary  is  first  recorded  by  a  set  number  cf  talkers,  usually  5  or 
10.  A  total  of  thirty  different  preordained  word  orders  are  available 
to  eliminate  possible  learning  effects  by  listening  subjects.  The 
microphone  used  to  make  the  recordings  and  the  ambient  acoustic 
environment  are  selected  according  to  the  proposed  operational 
deployment  of  the  communication  system  under  test.  The  recordings  are 
then  processed  by  passing  them  through  the  actual  communication  system 
(or  a  laboratory  simulation  if  this  is  not  possible)  and  re-recorded. 
The  processed  word  lists  are  then  presented  aurally  to  a  panel  of 
listeners  via  the  appropriate  earphone  transducer  which  is  usually 
contained  in  a  helmet,  headset  or  handset.  At  the  same  time  the 
listeners  are  also  presented  the  appropriate  word  pair  visually  on  a 
VDU.  They  then  select  the  word  they  thought  they  heard  by  pressing  one 
of  two  buttons.  The  intelligibility  score  is  expressed  as  a  proportion 
of  the  number  of  words  correctly  identified  by  the  listeners,  with  the 
scores  adjusted  for  chance  so  that  a  subject  guessing  will  score  0% 
rather  than  50%. 

To  ensure  adequate  stability  of  the  results,  tests  are  conducted  using 
not  less  than  five  talkers  and  eight  listeners.  Each  of  the  talkers 
utters  the  192  word  vocabulary  at  a  rate  of  one  word  every  1.33 
seconds;  a  five  talker  test  therefore  lasts  approximately  twenty-five 
minutes.  Male  and  female  subjects  are  recruited  from  the  local 
population  and  paid  for  their  participation  as  listeners  in  the  tests 
which  are  normally  conducted  three  mornings  a  week.  The  implcmontstion 
is  very  close  to  that  employed  by  US  testing  agencies,  thereby  allowing 
the  exchange  of  speech  material  between  the  US  and  the  UK  for 
comparative  assessment. 


1.2  RSRE  Implementation 

The  DRT  system  developed  at  RSRS  involves  a  network  of  British 
Broadcasting  Corporation  (BBC)  microcomputers  comprising  12  listener 
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stations,  a  file-server  and  a  master  station,  the  listeners'  VDUs  being 
installed  in  two  high-noise  chambers.  The  listeners  indicate  their 
choice  of  word  by  pressing  either  the  left  or  right  button  on  a  box 
connected  to  their  microcomputer  and  the  responses  are  logged  at  the 
master  station.  At  the  end  of  a  listening  session  the  results  for  all 
the  listening  subjects  are  amalgamated  and  stored  in  a  text  file  known 
as  a  "condition"  file.  This  file  contains  header  infoxrmation 
describing  the  test  environment  and  blocks  of  user  responses;  each 
being  signified  by  a  *1'  for  a  correct  response  and  a  'O'  for  an 
incorrect  response. 

The  results  can  be  analysed  locally  on  a  BBC  microcomputer  by  running 
the  program  DRTSTAT,  however  the  full  package  is  available  on  a-DEC  VAX 
hardware  in  N  building  South  Site.  The  terminal  in  the  Acoustic 
Laiboratory  connected  to  the  VAX  is  a  BBC  microcomputer  running  terminal 
emulation  software.  This  means  that  the  results  held  on  a  BBC 
formatted  disc  can  be  sent  down  the  communications  link  to  create  text 
files  on  the  VAX,  which  then  act  as  input  to  the  Analysis  of  Variance 
prograun,  DRTSTAT. 


1 . 3  Statistical  Analysis 

The  DRT  statistics  package  consists  of  three  programs  namely  DRTSTAT, 
DRTFILE  and  DRTEDIT,  which  are  run  by  typing  SDRTl,  @DRT2  and  @DRT3 
respectively  at  the  VAX  terminal. 

@DRT1  takes  a  previously  stored  condition  file  and  produces  statistical 
analysis  of  the  condition.  It  also  stores  the  basic  statistical 
(listener  and  talker  identification,  overall  mark,  attribute  and 
individual  listener  scores)  from  every  DRT  experiment  and  stores  uhe 
results  in  the  file  CONDMEANS.DAT. 

@DRT2  reads  CONDMEANS.DAT  and  takes  user  selected  conditions  and 
results  from  the  analysis  carried  out  by  6DRT1  and  reformats  them  for 
subsequent  analysis  by  the  program  PANOVUN.  This  Analysis  of  Variance 
computer  program  is  part  of  the  suite  of  statistical  analysis  software 
supplied  by  the  Institute  of  Aviation  Medicine  at  Farnborough  and  is 
available  on  the  VAX. 

@DRT3  provides  a  facility  to  list  and  edit  the  running  file  of  basic 
statistics  updated  by  @DRT2<, 


1.4  Summary  of  Output  from  DRTSTAT 

The  statistical  information  output  from  the  program  was  tailored  to 
suit  the  particular  requirements  of  the  DRT  experiments.  This  section 
describes  the  rationale  behind  the  selection  of  the  output.  See  Annex 
A  for  a  sample  printout. 
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1.4.1  Basic  Statistics 

The  fina]  result  to  be  output  is  the  overall  mean  score  for  the 
particular  condition  being  tested  along  with  the  standard  error  of 
listeners  scores  when  meaned  over  talkers.  However  it  is  also  useful 
to  be  able  to  see,  in  tabular  form,  the  performance  of  each  individual 
listener/talker  combination.  This  is  output  under  the  heading 
Individual  Scores. 

As  has  already  been  mentioned,  the  DRT  tests  the  condition  for  six 
major  attributes  of  speech.  An  output  is  therefore  required  which  will 
show  condition  performance  when  each  attribute  was  present  and  absent 
in  the  utterance.  This  is  output  under  the  heading  Attribute  Scores. 


1.4.2  Analysis  of  Variance 

Two  way  Analysis  of  Variance  is  designed  to  show  the  main  and 
interactive  effects  of  two  factors,  each  factor  having  more  than  one 
level.  In  the  case  of  DRT  it  was  decided  that  the  two  main  factors 
were  Listeners  and  Talkers,  with  the  list  parts  one  and  two  being  the 
replications.  If  the  talker/listener  interaction  is  shown  not  to  be  a 
significant  effect  then  th  Newman-Keuls  test  can  be  applied  to 
listeners  and  talkers  separately. 


1.4.3  Tukev  Test  for  Nonaddativitv 

If  the  experiments  were  to  be  carried  out  on  a  single  list  (i.e. 
without  replicaticn) ,  then  the  Tukey  test  for  nonaddativity  would  be 
used  to  test  the  presence  of  a  Talker  X  Listener  interaction.  If  the 
test  does  not  indicate  the  presence  of  an  interaction,  then  the 
analysis  of  the  main  effects  is  carried  out  as  described  in  the 
preceding  section. 
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2 

Note: 


OPERATING  INSTRUCTIONS 


@DRT1  requires  the  printer  to  be  put  on-line  for  hard  copy 
before  it  is  run. 

@DRT2  will  display  messages  as  the  required  files  are  generated 
and  allow  immediate  running  of  PANCVUN  if  desired. 


@DRT1  Enter  Condition  Filename 


@DRT2  Enter  Control  Filename 


Enter  the  name  of  the 
condition  for  analysis  in 
the  form  GND9999 

Enter  the  neune  of  the  file 
containing  the  control 
information.  (See  Annex  A. 4 
for  example) 


@DRT3  1  Display  Index  Select  option  required 

2  Delete  an  Entry 

3  End 
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3  PRCXJRAM  SPECIFICATIONS 
3.1  0DRT1  Prooram  Specification 

It  was  decided  that  a  total  of  five  ^etB  of  analysis  be  carried  out  on 
the  results  of  the  DRT  tests  and  their  implementation  along  with  a 
basic  description  is  described  below. 


3.1.1  Tukev  Test  for  Nonaddativitv 

If  there  is  only  one  observation  in  each  cell  of  a  p*q  factorial 
experiment  there  can  be  no  within-cell  variation  and  hence  no  direct 
estimate  of  experimental  error.  The  following  two  models  may  be 
postulated  to  underlie  the  observed  data: 


(1)  Xjj  =  p  +  a,  +  Pj  +  £(, 

(2)  Xy  =  p  +  a,  -s-  Pj  +  aPij  +  Cy 

In  (1)  no  interaction  effect  is  postulated;  hence  all  sources  of 
variation  other  than  main  effects  are  considered  to  be  part  of  the 
experimental  error. 

In  (2)  an  interaction  term  is  postulated.  This  interaction  term  may  be 
considered  as  a  measure  of  nonaddativity  of  the  main  effects. 

Tukey  developed  a  test  applicable  to  the  case  in  which  there  is  a 
single  observation  per  cell.  The  test  is  called  a  test  for 
'nonadditivity'  and  its  purpose  is  to  help  in  the  decision  between 
models  (1)  and  (2).  In  his  approach  to  the  problem,  Tukey  starts  with 
the  model: 


(3)  Xy  «=  p  +  a,  +  P,  +  Xa,  P,  +  Cy 


where  a,  Pj  is  the  product  of  the  main  effects  and  1  is  a  regression 
coefficient.  here  the  product  term  is  that  part  of  the  interaction 
which  can  be  expressed  as  the  product  of  main  effects.  A  test  in  the 
hypothesis  1  ■>  0  is  equivalent  to  a  test  on  the  hypothesis  that  the 
product  terms  of  this  form  do  not  contribute  to  the  prediction  of  the 
Xy. 

In  applying  Tukey 's  test  for  nonadditivity  to  a  p*q  factorial 
experiment  with  one  observation  per  cell,  the  analysis  of  variance 
takes  the  form: 
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Source  of  Variation 

SS 

d.f. 

MS 

A 

SSa 

p-1 

HSa 

B 

SSb 

q-1 

MSb 

Residual 

SSres 

(p-l)(q-l) 

Nonadditivity 

SSnonadd 

1 

MSnonadd 

Balance 

SSbal 

(p-l)(q-l)-l 

MSbal 

where  SSbal  SSres  -  SSnonadd 
and  SSrea  *  SStotal  SSa  -  SSb 


The  test  for  nonadditivity  is  given  by: 

F  =  MSnonadd  /  MSbal 

When  this  F  ratio  exceeds  the  corresponding  critical  value,  the 
hypothesis  that  model  (1)  (ie  X  =  0)  is  appropriate  is  rejected, 
Tukey's  test  for  nonadditivity  is  sensitive  only  to  one  source  of 
nonadditivity  -  that  associated  with  a  component  of  interaction 
represented  by  a,  Pj. 

For  full  description  of  test  see  Winer  -  'statistical  Principles  In 
Experimental  Design'  pp  394-397. 

Assuming  a  matrix  of  n  listeners  (i)  by  m  talkers  (j)  and  results  X,), 
the  equations  are  as  follows: 


SSa  -  (1/m)  *  ( sum  i  )*  -  (l/(n*m) 

(-•I 

■ 

SSb  •-  (l/n) '  (  sum  j  -  (l/(n*CT) 


SStotal^Y,  E  -  (!/(«*«)  )*  (E  E 

I  J  \  I  J 


SSres  -  SSab  “  SStotal  -  SSa  -  SSb 


r(EEx,f 


SSnonadd 


((Ec.».r)/((Ec;).(Ec;)) 


SSbal  •»  SSres  -  SSnonadd 
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MSx  =  (sax)  /  (d.f.x) 


Computational  symbols 
( 


(6) 

=  vc  E  E 

1-1  w-l  ) 

(7) 

=  1/r  E  (  E 

y-i  \i-i  / 

(8) 

-  l/(r*c)  (EE^f^ 

\  1  J 

r 

(9)  =EE^^ 

1  ; 

(10) 

«  1/c  E 

(10) 

-  l/(r*c)  E  E 

1  J 

(12) 

-  l/r  E 

1- 1 

(13) 

“  E 

J 

Then: 

SSlisners  •=  (6) 

(8) 

SStalkers  =  ( 7) 

(8) 

SStotal  -  (:.-) 

(8) 

SSres  «•  SStotal  -  SSlisners  -  SStalkers 


3.1.2  Two  Wav  Analysis  of  Variance 

This  analysis  allows  us  to  examine  the  combined  effect  of  the  two 
factors  'talkers'  and  'listeners'  ie  the  interactive  effect. 
Replications  in  this  case  refer  to  parts  1  and  2  of  the  talker  list. 
The  form  of  the  Analysis  of  variance  is  as  follows: 


Source 

SS 

d.f. 

MS 

F 

Rows 

SSr 

(r-1) 

SSr/ (r-1) 

HSr/MSerror 

Columns 

SSc 

(c-1) 

SSc/ (c-1) 

MSc/MSerror 

RxC 

SSrc 

(r-l)(c-l) 

SSrc/(r-l) (c-1) 

MSrc/MSerror 

Error 

SSerr 

rc(n-l) 

SSerr/rc(n-l) 

Total 

SStotal 

rcn-1 

Computation  of  the  various  values  is  as  follows: 
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c.^itt 

•'  I  j  k 


G^ttt 
I  J  » 


Computational  symbols 

(1) 

- 

cVrcn 

(2) 

-  tttxl 

t  J  k 

(3) 

ts 

(  E  R? )  /  cn 

t 

(4) 

-  (^c])/rn 

(5) 

EE  1  EXyj  j  /n 

ie  square  each  cells 

total  and  sum 

Then; 

SSr 

a 

(3) 

-  (1) 

SSc 

a 

(4) 

-  (1) 

SSrc 

a 

(5) 

-  (3)  -  (4)  +  (1) 

SSerror 

a 

(2) 

-  (5) 

SStotal 

a 

(2) 

-  (1) 

3.1.3  Nevman-KeulB  test 

The  Nevnnan-Xeuls  test  for  the  comparison  of  treatment  means  uses  the 
■Studentised  Range  Statistic',  Qr.  where  r  is  the  number  of  steps 
separating  two  moans  on  an  ordered  scale. 

Further  explanation  of  the  test  will  be  given  in  terms  of  the  DRT 
experiment  involving  six  talkers  and  weight  listeners.  We  shall  for 
the  moment  consider  the  talkers,  and  investigate  the  difference  between 
their  means. 


8 


posn  a  b  c  d  e  f 


Fig  3.1 


The  Talker  meano  ohould  first  be_ranked  T{1)~*T{6)  where  T{1)  is  the 
smallest  and  T(6)  the  largest.  T(6)  is  defined  as  being  six  steps 
from  T(l)  .  In  general  the  number  of  steps  between  the  ordered  means 
is  j-i+1. 

Thus,  06  -  [  r(6)  -  r(l)  )  /  'iaserror/n 


n  -  8  here  (16  if  2  reps) 

There  is  only  one  Q6  statistic.  There  are  however  two  Q5  statistics, 
neunely: 


05  -  (  T(6)  -  ?(2)  )  /  ^NSerror/n 


and 


05  -  I  r(5)  -?(!))  /  ^HSerror/n 


Critical  values  for  Qr  are  obtained  from  teddies  of  the  studentised 
range  statistic  by  setting  r  equal  to  the  range.  For  example,  the 
critical  value  for  Q6  in  a  0.01-level  test  is  Q.99  (6,f)  where  f  is  the 
degrees  of  freedom  for  MSerror. 

In  making  a  large  number  of  test  it  is  more  convenient  to  obtain  a 
critical  value  for  the  difference  between  two  means: 


T{6)  -  T(l)  --QR  yj  MSerror  In 

The  six  talker  means  are  designated  by  the  symbols  'a*  through  to  ‘f, 
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and  a  new  table  is  constructed  with  the  talker  means  in  ranked  order 
across  the  top. 


_ 

_ 

2; 

T3 

T2 

T3 

T3 

T3 

T3 

•  say ' 

c 

a 

d 

b 

f 

e 

r 

val 

c 

- 

Ti  -r3 

T4-T3 

T2  -T3 

re  -r3 

T5  -T3 

6 

xl  ,  say 

a 

- 

T^-Tl 

T2  -T1 

TS  -ri 

T5-T1 

5 

x2 

d 

- 

T2  -T4 

T6  -T4 

T5-T4 

4 

b 

- 

T6-T2 

T5  -T2 

3 

f 

- 

T5  -T6 

2 

where  val  Q.99  (r,f)  ^  MSerror/n 


Fig  3.2 


In  general  the  entry  in  column  j,  row  i  is  Tj-T,  .  Only  the  entries 
shown  in  this  table  need  be  computed.  It  should  be  noted  that  entries 
on  the  same  diagonal  on  the  table  have  the  seune  number  of  steps  between 
each  pair  of  means  being  compared. 

There  is  a  prescribed  sequence  in  which  tests  on  differences  between 
ordered  means  should  be  made: 

1  The  first  test  is  made  on  the  difference  in  the  upper  right- 

hand  corner.  This  is  the  difference  for  which  r  has  its 
maximum  value,  which  is  in  the  case  r  »  6.  Let  the  critical 
value  for  a  0.01-level  test  be  xl.  If  T5-TI  , 

ie  r(6)  -T^l)  does  not  exceed  this  value  one  need  not  carry  out 
any  further  tests  since,  if  there  is  no  significant  difference 
between  the  means  furthest  apart  from  each  other,  there  will 
clearly  be  no  significant  difference  between  any  pair  of  means 
closer  together. 

2  Tests  are  now  made  on  all  differences  for  which  r  =  5.  If  any 
of  the  difference  values  on  the  diagonal  for  which  r  »  5  do  not 
exceed  the  critical  value,  say  x2,  no  further  tests  are  made  on 
the  triangular  region  defined  by  the  rows  and  columns  of  which 
this  particular  difference  forme  the  upper  right-hand  corner. 
For  example,  if  TS-Tl  <  x2  there  would  be  no  need  to  test  the 
following: 
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M  -n 

T2-T1 

re-Ti 

T2  -TA 

Te-TA 

TB-TA 

T6-T2 

T5-T2 

T5-T6 

other  entries  on  the  diagonal,  though,  should  be  tested, 
eg  T6  -T3  ,  and  the  same  eliminating  procedure  carried  out 
where  necessary.  Thus  ws  work  through  on  the  entries  on  the 
table  given  above. 

The  schematic  representation  of  the  significant  differences  between 
means  of  talkers  take  the  form. 

c  a  d  b  f  e 


Fig  3.3 


The  letters  underlined  by  a  common  line  are  not  significantly  different 
to  each  other,  ie  there  is  no  significant  difference  between  15  ,  TS  , 
T2  and  TA  ,  nor  between  T3  ,  T1  ,  TA  . 


3.1.4  Attribute  Scores 

The  DRT  experiments  are  designed  to  examine  the  six  major  attributes 
which  apply  to  consonant  sounds  of  speech.  These  are: 


Voicing 

Nasality 

Sustention 

Sibilation 

Gravoness 

Compactness 


A  table  of  mean  and  standard  error  of  listeners  scores  is  therefore 
required  in  the  case  where  each  of  the  six  atcribotes  were  separately 
present  and  absent.  Also  included  in  the  tab7.e  is  the  mean  and 
standard  error  of  the  combined  present  and  absent  scores. 


3.1.5  Individual  Scores  and  Overall  DRT  Score 

The  basic  statistics  of  individual  listeners  scores  are  required  for 
detailed  examination  of  listeners  performance  and  to  allow 
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interpretation  of,  eg  the  Newman-Keula  test.  The  results  are  looked  at 
from  two  'directions',  ie  listener's  score  over  talkers,  and  talker 
score  over  listeners. 

The  final  statistic  required  is  the  overall  DRT  score  and  standard 
error. 

To  allow  for  archiving  of  results  for  further  analysis  certain  values 
are  stored  in  an  indexed  file,  the  key  field  being  the  condition 
number.  The  values  stores  are  as  follows: 


Condition  name 
Listener's  ID 
Listener's  mean  scores 
Listener's  attribute  scores 
Talker  results 
Talker's  ID 
Number  of  talkers 


3.2  0DRT2  Program  Specification 

In  order  to  utilise  the  results  saved  by  @DRTI,  a  program  is  required 
which  will  manipulate  selected  results  into  a  form  which  can  be 
processed  by  the  statistics  program  PANOVUN.  This  format,  known  as  EDA, 
is  an  Institute  of  Aviation  Medicine  standard. 


Inputs: 

1  Name  of  EDA  file  to  be  created. 

2  Attributes  required  from  conditions. 

3  Conditions  to  be  considered. 

4  Particular  listeners  of  interest. 

5  Format  of  EDA  file. 

6  Design  statements  required  for  PANOWN 

7  List  of  all  possible  listener  idents. 

Outputs: 

1  Data  output  files. 

1  Messages  to  screen. 
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3.3 


@DRT3  Program  Specification 


As  the  basic  statistics  from  @DRT3  are  automatically  stores  in  a 
running  file,  it  is  necessary  to  be  able  to  examine  that  file  to  see 
what  Condition  records  are  held.  It  is  also  necessary  to  be  edale  to 
delete  records  in  order  to  stop  the  file  from  becoming  too  large  and 
overloaded  with  outdated  information. 

A  program  is  required  which  will  satisfy  these  requirements  and  also  be 
written  in  such  a  way  that  will  allow  further  editing  facilities  to  be 
added  in  the  light  of  changes  in  policy. 


Input: 

Either  select  Index  of  Conditions  or  Condition  for 
deletion. 


Output: 

Either  list  of  Condition  numbers  or  notification  of 
Condition  deleted. 
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A1 


PROGRAM  DESCRIPTION 


Al. 1  gPRTl  Program  Description 

The  command  file  prompts  for  the  Condition  file,  copies  it  to  DRTRES 
and  runs  the  program  DRTSTAT.  As  can  be  seen  by  reference  to  the 
structure  chart,  the  program  consists  of  eleven  procedures  which  are 
described  as  below: 

INITIALISE 

1  Set  'Present  or  Absent’  array  to  zero. 

2  Set  listener  means  array  to  zero. 

3  Set  listener  ID  array  to  ' - •  and  talker  ID  array  to  ' - ’. 

UPDATECONDMEANS 

1  Collect  condition  name,  listener  and  talker  IDs,  attribute 

scores,  overall  means,  and  talker  results. 

2  Open  file. 

3  If  not  stored  by  a  previous  run,  store  results  using  condition 

as  key. 

4  Close  file. 

READSEQFILE 

1  Each  DRT  list  has  a  differer.  sequence  in  which  the  particular 

attribute  being  tested  is  present  or  ibsent.  These  sequences  are  read 
in  from  file  into  array  SEQ. 

READNKCONS 

1  Read  the  file  of  constants  required  by  the  Newman-Keuls  test. 

READWRITE 

1  Read  the  number  of  lines  in  parameter  from  DRTRES  and  output  to 

printer. 

READRESFIL 

1  Reset  DRTRES  for  reading 

2  Read  experiment  ID 

3  Read  version  number  and  output  necessary  number  of  lines. 

4  Build  up  array  of  talker  IDs. 
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5  If  DYNASEQ  was  used  then  for  each  of  the  presentations  in  end 
file,  calculate  which  SEQ  applies. 

6  Build  up  array  of  listener  IDs. 

7  Now  that  the  first  page  of  information  has  been  read,  the 
blocks  of  listener  response  can  be  read.  The  nesting  is  listeners 
within  parts  within  talkers.  Two  results  are  calculated;  the  attribute 
scores  and  the  overall  scores.  Attributes  are  only  calculated  if 
DYNASEQ  is  used.  In  which  case,  if  a  correct  response  is  read  (ie  a 
*1')  the  sequence  array  is  checked  to  find  if  the  attribute  was  present 
or  absent  and  the  correct  element  in  PORA  incremented.  The  variable 
ATT  cycles  through  the  six  attributes.  Both  the  attributes  and  the 
overall  results  have  the  algorithm  applied  which  gives  a  score  of  zero 
for  completely  random  responses. 

SINGLEREPS 

1  Produces  a  matrix  with  the  parte  average  for  the  Tukey  test. 

SETUNLNS 

1  Build  up  the  array  ULNS  which  gives  the  number  of  dashes 

required  for  up  to  12  positions  from  a  left  margin,  ie  pos  1  requires 
1  dash,  pos  3  requires  5.  This  is  required  when  deciding  how  many 
elements  in  the  Newman-Keuls  test  have  to  be  underlined. 

NKCALCS 

1  Read  in  constants  for  calculating  critical  values. 

2  Build  up  the  array  POSITIONARR. 

3  If  looking  at  talkers,  CONV  will  contain  numoflisners  and  vice- 
versa. 

4  Sort  the  array  of  means  (either  listener  or  talker)  into  lowest 
to  highest.  Also  sort  the  corresponding  position  in  POSITIONARR. 

0  Build  up  the  values  corresponding  to  Fig  3.2. 

6  If  no  significant  difference  as  in  3.3.1  the  print  out  ordered 
means  else  . . . 

7  Carry  out  comparisons  and  calculate  line  lengths  rec[uired 
(Fig  3.3). 

8  Printout. 


STATSONATTS 

1  Calculate  attribute  means,  sum  of  squares,  and  standard  error. 
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2  Write  results  for  attributes.  j 

ANOVAR 

1  Call  READRESFIL;  SINGLEREPS; 

2  Calculate  components  required  as  in  3.1  and  3.2. 

3  Write  results  for  Tukey  and  ANOVA, 

4  Calculate  listener  and  talker  means  and  call  NKCALCS. 

f 

5  If  DYNASEQ  was  used  then  call  STATSONATTS.  i' 

i 

6  Calculate  individual  means,  grandmean,  standard  deviation,  | 

standard  error  for  both  listeners  over  talkers  and  talkers  over 
listeners. 

7  Write  results  for  basic  stats. 


A1.2  gDRT2  Program  Description 

The  command  file  prompts  for  the  control  file  to  be  used  and  the  name 
of  the  DRTFILE  output  file  to  be  used  as  input  to  PANOVUN.  It  then 
runs  DRTFILE  followed  by  PANOVUN  if  desired.  Whilst  running,  DRTFILE 
will  report  any  errors  it  sees  in  the  input  data  but  whenever  possible 
will  allow  the  output  file  to  be  created.  This  is  to  allow  the 
operator  to  ignore  apparent  errors  for  a  particular  requirement.  As 
can  be  seen  on  the  structure  chart,  DRTFILE  consists  of  nine  procedures 
and  two  functions  as  follows: 

Function  INTTOCHR 

Takes  integers  in  the  range  0-99  and  turns  them  into  strings 
two  characters  long. 

Function  CHRTOINT 

Takes  strings  of  two  characters  and  returns  integers  in  the 
range  0-99. 

Function  TOCHARS 

This  function  simply  turns  the  real  numbers  stores  on  the  file 
CONDHEANS  into  characters  suitably  formatted  for  the  output 
text  file. 

1  If  a  score  of  100%  is  encountered  then  write  'lOO'elre.. 

2  Check  for  sign  to  be  inserted  in  first  element  of  field. 

3  Split  the  number  into  its  character  components. 
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Pack  in  into  a  string. 


EKTERNAMES 

1  Read  the  listener  ids  at  the  start  if  each  data  block 

and  check  to  ensure  that  each  4  character  name  exists  in 
the  set  of  listener  ids  at  the  end  of  the  control  file. 
Build  up  an  index  of  the  listener  id  numbers. 


UPDATE 

Using  recursion,  update  the  parameters  with  the  nested  values 
for  output. 

OUTPTNESTING 

Call  UPDATE  to  find  the  latest  values  for  the  nestings  and 
output  them. 

TOEDAFORM 

Read  the  control  file  to  find  the  format  required  for  the  EDA 
type  file  required  by  PANOVUN.  After  checking  for  errors  and 
extracting  the  required  information,  write  the  header 
information  to  the  output  file.  Cycle  through  the  temporary 
file  created  earlier  reading  the  listener  names,  extracting  the 
required  listeners  results,  checking  for  completeness,  sorting 
and  printing  out  the  results  to  the  output  file.  Finally  add 
the  required  number  of  -I's. 

DOEDAINPT 

Read  that  part  of  the  control  file  which  contains  the  design 
statements  for  PANOVUN  and  write  them  to  DRTEDA. 

OUTDATA 

1  Unpack  string  and  pick  our  3rd  component. 

2  Apply  a  hashing  function  to  obtain  a  corresponding 
number  for  the  attribute. 

3  Relate  the  has  number  to  the  required  value  in  PAORM  and 
write  the  value  after  applying  TOCHARS. 

SELECTDATA 

1  Krite  condition  number  and  listener  IDS. 

2  Step  through  the  string  of  attributes  required  and  call 
OUTDATA  for  each  until  complete. 
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RETRIEVEREC 


1  Locate  required  condition  using  end  number  as  key. 

2  If  data  located,  generate  attribute  means  and  call 
SELECTDATA  to  pick  out  required  elements. 

DECODESTRGS 

1  Find  length  of  user  input  strings  for  conditions  and 
attributes  required. 

2  Pick  out  each  condition  number  and.. 

3  If  a  range  has  been  specified,  work  out  what  each  number 
would  be  in  the  range  and  call  RETIEVEREC  for  each  else.. 

4  Call  RETRIEVEREC  for  condition  number. 

Main  program 

1  Open  files. 

2  Read  input  file  and  build  up  strings  of  attributes 
required  and  conditions. 

3  Issue  running  message. 

4  Call  OECODESTRINGS. 

5  If  any  of  the  conditions  required  do  not  exist  then 
report  and  stop  olse... 

6  Call  TOEDAFOPM  to  take  DRTTHPFLl  and  reformat  it  to 
DRTTMPFL2 . 

7  If  any  missing  values  exist  then  report. 

8  Generate  file  of  statements  required  for  PAMOVUN. 

A1.3  PDRT3  Program  Descrlution 

The  conunand  file  runs  the  program  DRTEDIT.  This  program  only  requires 
three  procedures,  ie  DISPINDEX,  DELENTRY  and  DISPHENU. 

1  DISPINDEX  opens  the  file  for  sequential  access,  reads  each 
entry  and  displays  the  condition  number. 

2  DELENTRY  opens  the  file  for  keyed  access,  find  the  required 
entry  and  deletes  it  from  the  file. 

3  DISPHENU  displays  the  options  and  the  program  CASE  statement 
vectors  to  the  required  procedure. 


DRT1.COM 

$  WR:=WRITE  SYS$0UTPUT 
$  WR  "  " 

S  WR  "DRT  STATISTICS  PROGRAM" 

S  WR  " - - 

$  WR  "  " 

$  AGAIN: 

$  INQUIRE  FILNAM  "Enter  filename" 

$  SET  NOON 

$  COPY  'FILNAM'  DRTRES.DAT 
$  IF  .NOT.  $STATUS  THEN  GOTO  AGAIN 
$  SET  ON 

$  PURGE  DRTRES.DAT 
$  RUN  DRTSTAT 

DRT2.COM 
S  DEASSIGN/ALL 
$  SET  TERM/WIDTH=132 
S  SET  TERM/NOWRAP 
S  WR;*WRITE  SYS$OUTPUT 
S  WR  "  " 

$  WR  "DRT  RESULTS  REFORMATTER  PROGRAM" 

S  WR  " - " 

S  WR  "  " 

$  INQUIRE  FILNAM  "Enter  Control  Filename" 

$  ASSIGN  'FILNAM'  infile 
$  ASSIGN  DRTTHPFL2.DAT  outfile 
$  ASSIGN/USER^MODE  'F$LOGICAL( "TT" T  SYSSINPUT 
$  PURGE  DRTEDA.DAT 
$  PURGE  DRTTHPFL2.DAT 
S  PURGE  DRTTHPFL1.DAT 
$  RUN  DRTFILE 

S  INQUIRE  RUNPAN  "Run  Panovun?  (Y/N)" 

$  IF  RUNPAN  .EQS.  "N"  THEN  GOTO  done 
$  PURGE  DRTRESULTS.LIS 
S  ASSIGN  DRTRESULTS.LIS  FOR013 
$  DEASSIGN  outfile 
$  ASSIGN  'FSLOGICAL("TT") •  FOR014 
$  ASSIGN  DRTTHPFL2.DAT  FOR015 
S  ASSIGN  'FSLOGICAL("TT") •  FOROIS 
$  ASSIGN  DRTONELINE.DAT  FOR017 
$  ASSIGN  •F$LOGICAL("TT") '  FOR018 
S  ASSIGN  'FSLOGICAL("TT") •  FOR019 
$  RUN  PANOVUN 
$  TYPE  DRTRESULTS 
$  done: 

$  DEASSICN/ALL 
DRT3.COM 

$  ASSIGN/USER__MODE  SYSJCOMMAND  SYSSINPUT 
$  RUN  DRTEDIT 
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PROGRAM  DRTSTAT { INPUT , OUTPUT . DRTRES , DRTCONS , CONDMEANS , DYNASEQ ) ; 


(*  RSRE  Malvern  1985.  *) 
(*  This  program  takes  output  from  Diagnostic  Rhyme  Tests  *) 
(*  produced  by  Multi-Listener  Master  and  performs  the  *) 
(*  statistical  tests  defined  in  the  specification.  *) 
(*  Although  designed  initially  to  run  on  a  DEC  VAX  vers4,*) 
(*  it  has  few  Vax  dependent  features  to  allow  easy  *) 
(*  conversion  for  other  PASCAL  compilers.  *) 


CONST 

NUM0FREPS*2 ; 

TYPE 

STATABC4-ARRAY[1..12]  OF  REAL; 

ARR0FMEANS*ARRAY(1. .12]  OF  REAL; 

ARR0FUNLNS*ARRAY(1. .12)  OF  INTEGER; 

STR1NG=PACKED  ARRAY (1.. 7)  OF  CHAR; 

STROFFOUR=PACKED  ARRAY {1.. 4]  OF  CHAR; 

STROFTHREEaPACKED  ARRAY(1..3)  OF  CHAR; 

CONDRECaRECORD 

CONDITION: (KEY (0) )  PACKED  ARRAY (1.. 7)  OF  CHAR; 
LNAMS : ARRAY (1. .12)  OF  STROFFOUR; 

LI SNEROAMS ; ARROFMEANS ; 

ATTRIBUTES:ARRAY(1, .12,1. .2,1. .6)  OF  REAL; 
TLKRESMAT;ARRAY(1..12,1..12)  OF  REAL; 

TNAMS : ARRAY (1. .12)  OF  STROFTHREE; 

TLKNUH: INTEGER ; 

END; 


VAR 

PORA:ARRAY(1..12,1..2,1..6)  OF  REAL; 
NUMOFTALKERS , NUMOFLISNERS , NUMOFPRES : INTEGER; 
DRTRES , DRTCONS , DYNASEQ : TEXT ; 

DONE: BOOLEAN; 

0AM , GRANDMEAN , STANDEV , STANERR , OASTANERR ; RE AL ; 
SSROK, SSCOL, SSROWCOL, SSERR, SSTOTAL, MSERR: REAL; 
OPTION, I , J, K: INTEGER; 

RESMAT : ARRAY (1. .12,1. .12,1. .NUMOFREPS)  OF  REAL; 
SINREPMAT:ARRAY(1..12,1..12)  OF  REAL; 
SEQS:ARRAY(1. .30,1. .96,1. .2)  OF  CHAR; 

RET; CHAR; 

UNLNS:ARROFUNLNS; 

LMEANS , TMEANS : ARROFMEANS ; 
ATTMEAN:ARRAY[1..2,1..6]  OF  REAL; 
PAMEAN:ARRAY(1. .6)  OF  REAL; 

NKCONS:STATABC4; 

EXPID: STRING; 

CONDMEANS: FILE  OF  CONDREC; 

OASTANDEV:REAL; 

LNRIDiARRAY(1..12)  OF  STROFFOUR; 
LNRIDARR;ARRAY(1..4)  OF  CHAR; 

TLKID ; ARRAY (1.. 12)  OF  STROFTHREE; 
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TLKIDARR;ARRAy[l. .3J  OF  CHAR; 
NVAL:INTEGER; 

ATTSS;ARRAY[1. .2,1. .6)  OF  REAL; 
ATTSE;ARRAY[1..3,1..6)  OF  REAL; 
TITLSTR; PACKED  ARRAY [1.. 66]  OF  CHAR; 
SUBTITL; PACKED  ARRAY (1.. 11)  OF  CHAR; 
TITLPTR;INTEGER; 

MATTSS:ARRAYtl. .6]  OF  REAL; 
SEQNUM:ARRAY(1. .30)  OF  INTEGER; 


PROCEDURE  INITIALISE; 
(*  Bet  arrays  to  0  *) 


BEGIN 

FOR  I;=l  TO  12  DO 
BEGIN 

FOR  J:=l  TO  2  DO 

FOR  K;=l  TO  6  DO 

PORA(I,J,K) :=0; 

LMEANS(I);=0; 

LNRID(I):  =  ' - '; 

TLKID( I)  ;  = 

END; 

END; 

PROCEDURE  UPDATECONDHEANS ; 

(*  add  overall  moan,  lioners, attributes  to  file  *) 

VAR 

COND:CONDREC; 

BEGIN 

WITH  COND  DO 
BEGIN 

CONDITION ;-EXPID; 

FOR  I:»l  TO  12  DO 
BEGIN 

I.NAMS(I):=LNRID(I); 

FOR  J:-l  TO  2  DO 

FOR  K:«l  TO  6  DO 

ATTRIBUTES[I,J,K):»PORA(I,J,K); 
LISNEROAHS [ I ] : -LKEANS ( I ) ; 

TLKNUM ; -NUHOFTALKERS ; 

FOR  J;=l  TO  NUHOFTALKERS  DO 

TLKRESHAT ( I , J ) : “SINREPHAT ) I , J ) 

END; 

FOR  I:»=l  TO  NUHOFTALKERS  DO  TNAMS(I)  ;»TLKIDtI)  ; 
END; 

OPEN(CONDMEANS, 

HISTORY! -UNKNOWN, 

ORGANIZATION:-INDEXED, 

ACCESS_HETHOD! -KEYED) ; 

RESETK ( CONDHEANS , 0 ) ; 

FINDK ( CONDHEANS , 0 , COND . CONDITION ) ; 
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IF  UFB(CONDMEANS) 

THEN 

VmiTE ( CONDMEANS , COND ) 

ELSE 

WRITELN (' entry  in  CONDKEANS.DAT  already  exists’); 
CLOSE ( CONDMEANS ) ; 

END; 

PROCEDURE  READSEQFILE; 

(*  create  matrix  of  DYNA  sequences  for  attributes  *) 

VAR 

I, J,K: INTEGER; 

CH:CHAR; 

BEGIN 

RESET (DYNASEQ) ; 

FOR  I;*l  TO  30  DO 

FOR  K:=l  TO  2  DO 
BEGIN 

FOR  J:«=l  TO  96  DO 
BEGIN 

READ ( DYNASEQ, CH); 

SEQS(I,J,K):=CH; 

END; 

READLN( DYNASEQ) ; 

END; 

END; 

PROCEDURE  READNKCONS; 

(*  read  file  of  constants  for  newman-)teulB  test  *) 

BEGIN 

RESET { DRTCONS ) ; 

FOR  I:«2  TO  12  DO  READLN( DRTCONS, NKCONS( I )) ; 

END; 

PROCEDURE  READWRITE (LINES; INTEGER); 

{*  read  lines  from  drtres.dat  and  output  to  printer  *) 
VAR 

I : INTEGER; 

CH;CHAR; 

BEGIN 

FOR  I:-l  TO  LINES  DO 
BEGIN 

WHILE  NOT  EOLN(DRTRES)  DO 
BEGIN 

READ(DRTRES,CH); 

WRITE (CH); 

END; 

WRITELN; 

READLN(DRTRES); 

END; 

END; 

PROCEDURE  READRESFIL; 

(*  read  drtres  and  build  up  matrices  of  results  *) 
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VAR 

CH , VERS , AORB : CHAR ; 

RIGHTS , WRONGS , L , OPLINES : INTEGER ; 

NUMTLK,TOTPERATT;INTEGER; 

ATT, SEQJ: INTEGER; 

BEGIN 

FOR  I:»l  TO  3  DO  WRITELN; 

RESET (ORTRES) ; 

FOR  I:*l  TO  14  DO  BEGIN  READ(DRTRES,CH) ;  WRITE ( CH ); END; 

READLN ( DRTRES , EXPID ) ; 

WRITELN (EXPID) ; 

READLN ( DRTRES ) ; 

WRITELN; 

(*  output  experiment  details  *) 

FOR  I:«l  TO  31  DO  BEGIN  READ ( DRTRES , CH ) ;  WRITE(CH);END; 

READLN (DRTRES, VERS) ; 

WRITELN (VERS); 

IF  VERS='l'  THEN  BEGIN 

READWRITE(21); 

FOR  I:=l  TO  15  DO 
BEGIN 

RE AD ( DRTRES, CH ) ; 

WRITE (CH) ; 

END; 

END 

ELSE 

BEGIN 

READWRITE(18); 

FOR  I:»l  TO  6  DO 
BEGIN 

READ(DRTFES,CH); 

WRITE (CH); 

ElO; 

END; 

READLN ( DRTRES , NUMOFPRES ) ; 

WRITELN (NUHOFPRSS  s  3 ) ; 

READWRITE ( 1 ) ; 

(*  read  list  numbers  and  calculate  which  row 

In  matrix  SEQ  each  presentation  corresponds  to  *) 
NUMTLK:"0;  SEQNUH(1] :>'1000; (*  tested  later  for  attributes  *) 
FOR  I;«l  TO  NUMOFPRES  DO 
BEGIN 

IF  I  MOD  2  <>  0 

THEN  (•  build  up  talker  id  strings  *) 

BEGIN 

NUHTLK:-NUMTLK+1; 

FOR  Jj»l  TO  3  DO 
BEGIN 

READ ( DRTRES, TLKIDARRIJ) ) ; 

VJRITE(TLKIDARR(J]) 

END; 

PACK(TLKIDARR, 1,TLKID(NUHTLKJ ) ; 

FOR  J;«l  TO  13  DO  BEGIN  READ ( DRTRES, CH ); WRITE (CH)  END 
END 
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ELSE 

FOR  J;=a  TO  16  DO  BEGIN  READ (DRTRES, CH) ; WRITE (CH)  END; 
READ(DRTRES,CH); 

WRITE ( CH); 

IF  CH= • . ' 

THEN 

BEGIN 

READ ( DRTRES, CH); 

WRITE (CH) 

END 

ELSE 

BEGIN 

SEQNUM ( NUMTLK ] : «ORD (CH)-48; 

READ(DRTRES,CH); 

WRITE (CH); 

SEQNUM ( NUMTLK ] : ■10*SEQNUM [ NUMTLK  J  +ORD ( CH ) -48 ; 

READ ( DRTRES , AORB ) ; 

IF  AORB*'B'  THEN  SEQNUM(NUMTLK) :*SEQNUM( NUMTLK) +15 ; 
WRITE(AORB); 

FOR  J:=l  TO  2  DO  BEGIN  READ(DRTRES,CH) ; WRITE (CH) ; END; 
END; 

READLN( DRTRES) ; 

WRITELN; 

END; 

READLN( DRTRES); 

FOR  I;«l  TO  10  DO  BEGIN  READ{DRTRES, CH) ;  WRITE(CH) ;END; 

READLN (DRTRES, NUMOFLISNERS) ; 

WRITELN (NUHOFLISNERS ; 3 ) ; 

READWRITE(l); 

(*  build  up  listener  ID  strings  *) 

FOR  I;-l  TO  NUHOFLISNERS  DO 
BEGIN 

FOR  Js«l  TO  4  DO 
BEGIN 

READ ( DRTRES , LNRIDARR ( J ) ) ; 

WRITE(LNRIDARR{JJ) 

END; 

PACK (LNRIDARR, 1 , LNRID ( I ) ) ; 

WHILE  NOT  EOLN  (DRTRES)  DO  BEGIN 

RE AD ( DRTRES, CH); 

WRITE ( CH ) ; 

END; 

READLN (DRTRES) ; WRITELN; 

END; 

READLN (DRTRES) ; WRITELN; 

WHILE  NOT  EOLN  (DRTRES)  DO  READ(DRTRES,CH); 

READLN ( DRTRES ) ; 

NUHOFTALKERSfNUHOFPRES  DIV  2; 

(*  read  0  and  I's  for  each  pres  and  build  up  results  matrix  *) 
FOR  J:»l  TO  NUMOFTALKERS  DO 
FOR  K 1-1  TO  NUMOFREPS  DO 

FOR  I:=l  TO  NUHOFLISNERS  DO 
BEGIN 

FOR  L;al  TO  2  DO 


All 


) 


BEGIN 

WHILE  NOT  EOLN  (DRTRES)  DO  READ(DRTRES,CH) ; 
READLN(DRTRES); 

END; 

RIGHTS :*0;  WRONGS:=0; 

SEQJ:*0;  ATT:*-!); 

FOR  L:*l  TO  3  DO 
BEGIN 

WHILE  NOT  EOLN  (DRTRES)  DO 
BEGIN 


SEQJ;»SEQJ+1; 

XTT;*ATT+1; 

IF  ATT»7  THEN  ATT:*1; 

READ(DRTRES,CH); 

IF  CH='l'  THEN  BEGIN 

RIGHTS :«RIGHTS+1; 

IF  SEQNUH[  1)01000  THEN 
IFSEQS(SEQNUM[J),SEQJ,K)='0' 
THEN  P0RA[I,1,ATT) := 
P0RA(I,1,ATTJ  +  1 
ELSE  P0RA(I,2,ATT) 

P0RA(I,2,ATT)  +  1; 

END 

ELSE  WRONGS !=WR0NGS+1; 

END; 

READLN( DRTRES); 

END; 


RESMAT  ( I ,  J ,  K  ) :  »=  ( ( RIGHTS-WRONCS )  *  100 )  / 


END; 


(RIGHTS+WRONGS); 


(■*  apply  conversion  to  attributes  *) 

TOTPERATT : -NUMOFFRES  •  8; 

FOR  I;»  1  TO  NUHOFLISNERS  DO 
FOR  J;*l  TO  2  DO 

FOR  K;»l  TO  6  DO 

PORA(I, J,K] :a((2*PORA(I,J,Kl~TOTPERATT)*100) 

/TOTPERATT; 

END; 


PROCEDURE  SINGLEREPS; 

(*  average  repetitions  for  tukey  test  *) 

BEGIN 

FOR  I:»l  TO  NUHOFLISNERS  DO 

FOR  J:«l  TO  NUHOFTALKERS  DO 

SINREPHAT ( I , J ] : » ( reshat { I , J , 1 J +RESMAT ( I , J , 2 ) ) /2 ; 

END; 


PROCEDURE  SETUNLNS; 

(*  calculate  underlining  for  NK  test  •) 
VAR  STEPS; INTEGER; 

BEGIN 

STEPS ;=1; 

FOR  I;=l  TO  12  DO 
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BEGIN 

UNLNS[I]:=STEPS; 

STEPS :=STEPS+2; 

END; 

END; 

PROCEDURE  NKCALCS(TREATMEANS:ARROFMEANS; 

TKLN: INTEGER) ; 

(*  newman-keuls  tast  for  liatenera  and  taDcers  *) 

LABEL 

99; 

VAR 

CONV , DASHES , LINESUB , NUMTKLN , PREVIOUS , TEHP2 : INTEGER ; 
SIGDIFF, SORTED: BOOLEAN; 

TEMP , Q9  9MULT : REAL ; 

NKMAT:ARRAY(1..12,1..13]  OF  REAL; 

LINELEN: ARRAY (1. .40)  OF  INTEGER; 

POSITIONARR : ARRAY (1. .12)  OF  INTEGER; 

BEGIN 

READNKCONS; 

NUMTKLN ;«TKLN; 

FOR  I:»l  TO  NUMTKLN  DO  POSITIONARR) I ): «I ; 

CONV ; -NUMOFTALKERS+NUMOFLI SNERS-NUMTKLN ; 

{*  sort  *) 

FOR  I;-l  TO  NUMTKLN  DO 
BEGIN 

SORTED: -TRUE; 

FOR  J:«l  TO  (NUMTKLN-1)  DO 

IF  TREATMEANS(J)>-TREATMEANS(J+1)  THEN 
BEGIN 

TEMP : -TREATMEANS ( J ) ; 

TREATMEANS ( J ) : -TREATHEANS (J+1); 
TREATMEANS ( J+1 ) : -TEMP; 
TEMP2;-POSITIONARR(J) ; 
POSITIONARR) J ) : -POSITIONARR ) J+1 ) 
POSITIONARR ) J+1 ) : -TEMP2 ; 

SORTED: -FALSE; 

END; 

END; 

(*  create  matrix  *) 

FOR  I:-l  TO  NUMTKLN-1  DO 

FOR  J:»l  TO  NUMTKLN-I  DO 

NKMAT ) I , J+I ) : -TREATMEANS ) J+I ) -TREATMEANS ) I ) 
(*  calculate  constant  table  multiplier  *) 
Q99MULT:-SQRT(MSERR/(NUHOFREPS*CONV) ) ; 

(*  aminend  matrix  to  contain  999  for  **  in  Winer  *) 
IF  NKMAT) l,NUMTKLN)<(NKCONS (NUMTKLN) *Q99MULT) 

(*  i.e.  if  no  eigdiff  in  top  right  of  matrix  *) 

(*  then  print  1  -  numtkln  and  underline  *) 

THEN 

BEGIN 

WRITE)’  ■); 

FOR  I:«l  TO  NUMTKLN  DO  WRITE ( POSITIONARR)!) :2); 
WRITELN; 
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WRITEC  *); 

FOR  I;=l  TO  NUMTKLN*2-1  DO  WRITEC-'); 

END 

ELSE 

BEGIN 

NKMAT ( 1 , NUMTKLN ] : *9 9 9 ; 

FOR  I; *2  TO  NUMTKLN-1  DO 
FOR  Jj*0  TO  I-l  DO 

IF  NKMAT[I-J,NUMTKLN-J)  > 

( NKCONS [ NUMTKLN+ l-I ) *Q99MULT ) 

THEN  NKMATII-J,NUMTKLN-J];*999; 

(*  0et  up  underlines  *) 

FOR  I;«l  TO  NUMTKLN-1  DO 
BEGIN 

LINELENtI]:-!; 

FOR  J;»I+1  TO  NUMTKLN  DO 
BEGIN 

(*  count  non-aigdiffs  *) 

IF  NKMAT(I,J)  <>  999 
THEN  LINELEN[I) ;=LINELEN(I1+1 
ELSE  GOTO  99; 

END; 

99:  END; 

(*  printout  *) 

WRITELN; 

WRITEC  ')! 

FOR  I:-l  TO  NUMTKLN  DO  WRITE(POSITIONARR( I ) : 2 ) ; 
WRITELN; 

SETUNLNS; 

DASHES s=l; 

PREVIOUS :■!; 

FOR  I;«l  TO  NUMTKLN- 1  DO 
BEGIN 

IF  LINELEN(IJ  <>  1  THEN 
BEGIN 

FOR  J;»l  TO  DASHES  DO  WRITEC  '); 

IF  (2*DASHES)+UNLNSILINELEN{I)  JOPREVIOUS  THEN 
FOR  J:«l  TO  UNLNS(LINELEN[I1  J  DO  WRITEC-'); 
WRITELN; 

END; 

PREVIOUS:-(2*DASHES)+UNLNS(LINELEN[I) ] ; 

DASHES ;»DASHES+1 ; 

END; 

END; 

END; 

PROCEDURE  CALCLMEANS; 

VAR 

MEANSTOT:REAL; 

BEGIN 

(*  calc  meana  for  N.K.  *) 

FOR  I:'=l  TO  NUMOFLISNERS  DO 
BEGIN 

KEANSTOT;=0; 

FOR  J:=l  TO  NUMOFTALKERS  DO 
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MEANSTOT:=MEANSTOT+SINREPMAT{I, J] ; 

LMEANS [ I ] : =MEANSTOT/NUMOFTALKERS ; 

END; 

END; 

PROCEDURE  CALCTMEANS; 

VAR 

MEANSTOT:REAL; 

BEGIN 

(*  calc  tmeana  for  N.K.  *) 

FOR  J:*l  TO  NUMOFTALKERS  DO 
BEGIN 

MEANSTOT:*0; 

FOR  I;*l  TO  NUMOFLISNERS  DO 

MEANSTOT;«MEANSTOT+SINREPMAT(I, J] ; 

TMEANS  [  J ]  •. -MEANSTOT/NUMOFLISNERS  ; 

END; 

END; 

PROCEDURE  STATSONATTS; 

BEGIN 

{*  state  on  attributes  *) 

FOR  J:=l  TO  2  DO 

FOR  K:=l  TO  6  DO 
BEGIN 

ATTMEAN(J,K] :=0; 

ATTSS(J,K] J=0; 

MATTSS(K] :=0 
END; 

FOR  J:«l  TO  2  DO 

FOR  K;»l  TO  6  DO 
BEGIN 

FOR  I:=l  TO  NUMOFLISNERS  DO 
BEGIN 

ATTMEAN(J,K):«ATTMEAN(J,KJ+PORA(I, J,K]; 

ATTSS(J,K)  ;-ATTSS(J,K]+SQR(POBA(I,J,K)); 

END; 

ATTMEAN [ J , K ] : «ATTMEAN [ J , K ) /NUMOFLISNERS ; 

END; 

(*  moan  of  P  or  A  *) 

FOR  K:*l  TO  6  DO 

PAME AN ( K ) : - ( ATTME AN  1 1 , X ) +ATTHEAN [ 2 , K ] ) / 2 ; 

(•  SS  of  P  or  A  *) 

FOR  Kj-1  to  6  DO 

FOR  I;»l  TO  NUMOFLISNERS  DO 
FOR  J;»l  TO  2  DO 

MATTSS(K];«MATTSS(K]+SQR(PORAlI,J,K]); 

{*  SE  Of  P  or  A  *) 

FOR  I;*l  TO  2  DO 

FOR  J:<=1  TO  6  DO 
IF  NUMOFLlSNERSal 
THEN 

ATTSE[I,J] J*0 
ELSE 

ATTSS 1 1 ,  J ) :  -SQRT { { ATTSS  [  I ,  J )  -NUMOFLISNERS*SQR (ATTMEAN  [  I ,  J ] ) 
)/ (NUMOFLISNERS* (NUMOFLISNERS-l) ) ) ; 
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{*  SE  of  means  of  P  and  A  *) 

NVAL : =NUMOFLI SNERS  *2 ; 

IF  NUM0FLISNERS>1  THEN 
FOR  J:=l  TO  6  DO 

ATTSE [ 3 , J ] : =SQRT ( (MATTSS [ J ] -NVAL*SQR ( PAKEAN [ J ] ) ) 
/(NVAL*(NVAL-1))); 

(*  printout  *) 

TITLSTR:^' Voicing  Nasality  Sustention  Sibilation  Graveness 

Compactness' ; 

TITLPTR:*!; 

FOR  I;»l  TO  4  DO  WRITELN; 

WRITELN ( '  ATTRIBUTE  SCORES ' ) ; 


VmiTELNC  - •); 

WRITELN;  WRITELN; 

WRITELN ( '  Present  Absent  Mean ' ) ; 

WRITELN ( '  mean  SE  mean  SE  mean  SE ' ) ; 

WRITELN; 


FOR  I:=l  TO  6  DO 
BEGIN 

SUBTITL : =SUBSTR ( TITLSTR, TITLPTR, 11 ) ; 

WRITELN (SUBTITL, '  ' , ATTMEANJ 1, I ) : 8: 1, ATTSE { 1 , I ) : 6 : 1, 
ATTMEAN(2,I] :8;1,ATTSE(2,I):6;1,PAMEAN(I) :8:1, 
ATTSE(3,I] :6:1) ; 

TITLPTR ; =TITLPTR+11 ; 

END; 

END; 

PROCEDURE  ANOVAR; 

{*  main  statistics  calculation  and  output  routine  *) 

VAR 

G , TOT3 , TOT4 , TOTS , TOT6 , TOT7 , TOT8 , TOT9 ; REAL ; 
TOT3SUBl,TOT4SUBl,TOT6SUB,TOT7SUB:REAL; 

GRANDTOT , SDTOT : REAL ; 

NONADDQUOT , NONADDSUBl , NONADDSUB2 , NONADDIV : REAL ; 
SSNONADD, SSBAL, HSBAL, SSRES :REAL; 

MSROW, ROWRAT,HSCOL, COLRAT,MSROWCOL, ROWCOLRAT; REAL; 

LI SNDF , TALKDF , LNTKDF , ERRDF , TOTDF , BALDF : INTEGER; 
COHPNT ; ARRAY (1. .12)  OF  REAL; 

COMPNTIO: ARRAY (1.. 12)  OF  REAL; 
TSSROW,TSSCOL,TMSROW,THSCOL;REAL; 

COMPNT12: ARRAY [1.. 12)  OF  REAL; 

SDLN,SDTKs ARRAY [1.. 12)  OF  REAL;  BEGIN 
Gs-0;COMPNT(2) :“0;TOT3:»0;TOT4;»0;TOT5i«0; 

FOR  I; “I  TO  NUMOFLISNERS  DO 
BEGIN 

TOT3SUB1:*=0; 

FOR  J;=l  TO  NUMOFTALKERS  DO 
BEGIN 

T0T3SUB1:*=T0T3SUB1 

+RESHAT ( I , J, 1 ) +RESMAT { I , J , 2 ) ; 
TOT5;«TOT5+ 

{RESMAT(I,J,1)+RESMAT(1,J,2))**2; 
FOR  K:=l  TO  NUHOFREPS  DO 
BEGIN 
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G:-G+RESMAT[I, J,K] ; 

COMPNT ( 2 ) :  *=COMPNT  [  2  ]  + 

RESMATtI,J,K]**2; 

END; 

END; 

TOT3 : *TOT3+TOT3SUBl**2 ; 

END; 

(*  build  component*  required  in  epecification  *) 
COMPNT [ 1 ] ; »G**2/ (NUMOFLISNERS^NUMOFTALKERS 

♦NUMOFREPS) ; 

COMPNT ( 3 ] : *TOT3 / ( NUMOFTALKERS  *NUMOFREPS ) ; 

COMPNT  1 5 ] : »TOT5 /NUMOFREPS ; 

TOT4;*0; 

FOR  J:»l  TO  NUMOFTALKERS  DO 
BEGIN 

TOT4SUB1:*0; 

FOR  I:«l  TO  NUMOFLISNERS  DO 
BEGIN 

TOT4SUB1 : «TOT4SUBl+RESMAT[ I , J, 1 ] 

+RESMAT[I, J,2) ; 

END; 

TOT4 ; =TOT4+TOT4SUBl*  *2 ; 

END; 

COMPNT( 4 ] : =TOT4/ (NUMOFLISNERS*NUHOFREPS ) ; 

S3ROW:=COHPNTt3)-COMPNT{l) ; 

SSCOL : »COMPNT [ 4 ] -COMPNT { 1 ] ; 

SSROWCOL : -COMPNT [ 5 ) -COMPNT ( 3 ) -COMPNT ( 4 ) +COMPNT ( 1 ) ; 
SSERR : -COMPNT ( 2 ) -COMPNT ( 5 ] ; 

SSTOTALs -COMPNT ( 2 1 -COMPNT ( i ) ; 

LISNDF ; -NUMOFLISNERS-1 ; 

TALKDF  J -NUMOFTALKERS- 1 ; 

LNTKDFs»LISNDF*TALKDF; 

ERRDF: -NUMOFLISNERS*NUMOFTALKERS* (NUHOFREPS-1 ) ; 
TOTDF ; -NUMOFLISNERS • NUMOFTALKERS  *  NUMOFREPS - 1 ; 
MSROW:»SSROW/LISNDF; 

MSCOL:»SSCOL/TALKDF; 

MSROWCOL: -SSROWCOL/LNTKDF; 

HSERR; -SSERR/ERRDF; 

ROWRAT : -MSROW/MSERR; 

COLRAT: -MSCOL/MSERR; 

ROWCOLRAT; -MSROWCOL/MSERR; 

(*  Tukey  Teat  for  Nonadditivity  *) 

SINGLEREPS; 

TOT6:-0;  TOT8:»0;  TOT9:»0; 

FOR  I 1-1  TO  NUMOFLISNERS  DO 
BEGIN 

TOT6SUBI-0; 

FOR  J:«l  TO  NUMOFTALKERS  DO 
BEGIN 

TOT6SUBl-TOT6SUB+SIfmEPMAT[I, J) ; 
TOT8:»TOT8+SINREPMAT(I,J] ; 

TOT9 : »TOT9+SINREPMAT( I , J ] **2 ; 

END; 

TOT6 : -TOT6+TOT6SUB*  *2 ; 
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COMPNTIO  [  I  ] ;  >=TOT6SUB/NUHOFTALKERS; 

END; 

COMPNT [ 6 ] : =TOT6/NUMOFTALKERS ; 

COMPNT[8] ;=TOT8**2/(NUMOFLISNERS*NUMOFTALKERS) ; 

COMPNT(9]:*TOT9; 

COMPNT [ 11 ] : *TOT8/ ( NUMOFLISNERS*NUMOFTALKERS ) ; 

TOT7:=0; 

FOR  J:-l  TO  NUMOFTALKERS  DO 
BEGIN 

TOT7SUB;»0; 

FOR  I:»l  TO  NUMOFLISNERS  DO 
BEGIN 

TOT7SUB : -TOT7  SUB+SINREPMAT ( I , J ) ;  1 

END;  , 

TOT7:*TOT7+TOT7SUB**2;  ' 

COMPNT12 [ J ] : «TOT7SUB/NUHOFLISNERS ;  ( 

END; 

COMPNT ( 7 ]; =TOT7 /NUMOFLI SNERS ; 

TSSROW:=COMPNT(6]-COMPNT(8] ; 

TSSCOL : =COMPNT ( 7 ) -COMPNT ( 8 ] ; 

SSRES : =COMPNT ( 9 ) -COMPNT ( 6 ] -COMPNT { 7 ) +COMPNT ( 8 ) ; 

(*  Calc  Nonadditivity  *) 

NONADDQUOT:=0;  NONADDSUB2:=0; 

FOR  I:=l  TO  NUMOFLISNERS  DO 
BEGIN 

NONADDSUBls=0  ; 

FOR  J;=l  TO  NUMOFTALKERS  DO 
BEGIN 

NONADDSUBl ; =NONADDSUBl+ (COMPNT12 ( J ) -COMPNT ( 11 ) ) 

*SINREPMAT(I,J); 

END; 

NONADDSUB2 : *NONADDSUB2+ ( COMPNTIO [ I ) -COMPNT [ 1 1 ) ) * • 2 ; 

NONADDQUOT : -NONADDQUOT+NONADDSUBl 

*(COHPNT10(I]-COMPNT(11) ) ; 

END; 

NONADDQUOT ; »NONADDQUOT  **2; 

NONADDIV;»0; 

FOR  J:«=l  TO  NUMOFTALKERS  DO 
BEGIN 

NONADDI V : -NONADDI V+ ( COMPNT12 ( J ) -COMPNT ( 1 1 ) ) *  *2 ; 

END; 

IF  NONADDIV  -  0  THEN  NONADDIV:-  0.001; 

HONADDIVj»NONADDSUB2*NONADDIV; 

SSNONADD : -NONADDQUOT/NONADD I V ; 

SSBAL : -SSRES-SSNONADD ; 

BALDFi-LNTKDF-1; 

TMSROW; -TSSROW/LISNDF; 

TMSCOL : »TSSCOL/TALKDF ; 

MSBAL;»SSBAL/BALDF; 

FOR  I:«l  TO  12  DO  TOITELN('  '); 

WRITELNC  TUKEY  TEST  FOR  NONADDATIVITY * ) ; 

WRITELNC  - •); 

FOR  I:=l  TO  3  DO  WRITELN; 

WRITELN (• Source  Sum  of  Degrees  of  Mean  '); 
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Square' ) ; 


WRITELN ( '  Squares  Freedom 

WRITELNC  •); 

WRITELN ( • LISNERS  ' , TSSROW: 13 : 2>LISNDF: 10 , TMSROW: 13:2); 

WRITEI.N( 'TALKERS  ' ,TSSCOL: 13 :2,TALKDF: 10,THSCOL; 13 ; 2 ) ; 

WRITELN { 'RESIDUAL' ,SSRES;13:2,LNTKDF:10) ; 

WRITELN (• NONADD  ' ,SSNONADD:13:2, '  1' ,SSNONADD:13;2) ; 

WRITELN ( 'BALANCE  ' ,SSBAL:13;2,BALDF:10,MSBAL:13:2) ; 

FOR  I;»l  TO  5  DO  WRITELN; 

WRITELNC  ANALYSIS  OF  VARIANCE'); 

WRITELNC  - •); 

WRITELN  C  • ) ; WRITELN ( '  ' ) ; 

WRITELN (' Source  Sum  of  Degrees  of  Mean'/ 

•  F-Ratio'); 

WRITELN ('  Squares  Freedom  Square'); 

WRITELNC  '); 

WRITELN( 'LISNERS  ' ,SSROW;13;2,LISNDF:10,MSROW:13:2,ROWRAT:11:2) ; 
WRITELN (' TALKERS  ' , SSCOL: 13 ; 2,TALKDF; 10,MSCOL: 13 s 2 , COLRAT: 11 : 2 ) ; 
WRITELNCLN  *  TK  '  , SSROWCOL:  13 ; 2,LNTKDF:  10,MSROWCOL:  13 : 2 , 

ROWCOLRAT:ll:2) ; 

WRITELN( 'ERROR  ' ,SSERR;13:2,ERRDF:10,MSERR;13:2) ; 

WRITELN ( 'TOTAL  ' ,SSTOTAL:13:2,TOTDF:10) ; 

(*  Newman  -  Keula  Teat  *) 

FOR  Is=l  TO  4  DO  WRITELN; 

WRITELN ( ’  NEWMAN  -  KEULS  TEST ' ) ; 

WRITELNC  - ');WRITELN; 

WRITELN ('LISTENERS  : ' ) ; WRITELN; 

NKCALCS ( LHEANS , NUMOFLISNERS ) ; 

WRITELN ; WRITELN; WRITELN ( 'TALKERS  : * ) ; 

NKCALCS ( TMEANS , NUMOFTALKERS ) ; 

IF  SEQNUH(1)<>1000  THEN  STATSONATTS; 

(*  calculate  treatment  means  for  basic  stats  *) 

SDTOT  s -0 ; GRANDTOT : -0 ; 

FOR  l;«l  TO  NUMOFLISNERS  DO 
BEGIN 

SDLN(IJ;»0; 

FOR  J:»l  TO  NUMOFTALKERS  DO 

SDLN(I] ;»SDLN(I)+SINREPMAT(I, J)**2; 

SDLN [ 1 ] : -SQRT ( ABS ( SDLN ( I J -NUMOFTALKERS  *(LHEANStI))**2) 

/ (NUMOFTALKERS-1 ) ) ; 

SDTOT  J -SDTOT+ { LMEANS ( I ) ) *  *2 ; 

GRANDTOT j  »GRANDTOT+LMEANS ( I ) ; 

END; 

GRANDMEAN  S -GRANDTOT/NUMOFLISNERS ; 

STANDEV: -SQRT ( ABS ( SDTOT-NUHOFLISNERS* GRANDMEAN*  *2 ) 

/  (NDMOFLISNERS-1)); 

OASTANDEV : -STANDEV ; 

STANERR : «STANDEV/S0RT ( NUMOFLISNERS ) ; 

FOR  I:«l  TO  10  DO  WRITELN; 


WRITELNC  INDIVIDUAL  SCORES'); 

WRITELNC  - •); 


FOR  I:»=l  TO  3  DO  WRITELN; 
WRITECLNR  Tl'); 

FOR  I; -2  TO  NUMOFTALKERS  DO  WRITE (' 


TM;1); 


WRITE { ’  SD ' ) ; 

WRITELN; 

FOR  I:=l  TO  NUMOFLISNERS  DO 
BEGIN 

WRITE(I:2); 

FOR  J;*l  TO  NUMOFTALKERS  DO  WRITE{SINREPMAT[ I, J] : 6: 1) ; 
WRITE(SDLN(I];6:1) ; 

WRITELN; 

END; 

WRITELN; 

WRITELN (’LNR  MEAN  OVER  TKR’); 

FOR  I:«l  TO  NUMOFLISNERS  DO 

WRITELN(I:2,LMEANS(I):13:1); 

WRITELN; 

WRITELN( 'STD.  DEVN.  «' ,STANDEV:S:1) ; 

WRITELN ( 'STD.  ERROR  =' ,STANERR:5:1) ; 

OASTANERR ; =STANERR ; 

SDTOT:=0; 

FOR  J:=l  TO  NUMOFTALKERS  DO 
BEGIN 

SDTK(J];=0; 

FOR  I:=  1  TO  NUMOFLISNERS  DO 

SDTK[J]  !-SDTK[J]+SINREPMAT(I,J)*'»2; 

SDTK(J) :-SQRT( (SDTK(J)-NUMOFLISNERS*TMEANS(J)»«2) 

/{NUMOFLISNERS-1)); 

SDTOT:-SDTOT+(TMEANS( J) )**2; 

END; 

STANDEV ; -SQRT ( ( SDTOT-NUMOFTALKERS*GRANDMEAN*  *  2 ) 

/  (NUMOFTALKERS-1) ) ; 

STANERR : -STANDEV/SORT ( NUMOFTALKERS ) ; 

FOR  I;-l  TO  5  DO  WRITELN; 

WRITE{'TKR  LI'); 

FOR  I;»2  TO  NUMOFLISNERS  DO 

IF  KIO  THEN  WRITEC  L',I;1)  ELSE  WRITE('  L',I:1); 
WRITE {'  SD'); 

WRITELN; 

FOR  I;»l  TO  NUMOFTALKERS  DO 
BEGIN 

WRITE(Ij2); 

FOR  J:»l  TO  NUMOFLISNERS  DO  WRITE{SINREPMAT(J,I]:6;1); 
WRITE{SDTK(I]:6:1); 

WRITELN; 

END; 

WRITELN; 

WRITELN ('TKR  MEAN  OVER  LNR'); 

FOR  Ifl  TO  NUMOFTALKERS  DO 

WRITELN(I:2,TMEANS(I];13;1); 

WRITELN; 

WRITELNCSTD.  DEVN.  «=',STANDEVj5;1); 

WRITELN (’STD.  ERROR  STANERR: 5:1) ; 

FOR  I:«=l  TO  5  DO  WRITELN; 

WRITELN (’  DRT  SCORE  »' ,GRANDMEAN;5:1) ; 

WRITELN (’STANDARD  ERROR  »' ,OASTANERR;5:l) ; 

FOR  I:«l  TO  5  DO  WRITELN; 
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END; 


PROCEDURE  BASIC_AND_ATT_STATS; 

BEGIN 

IF  SEQNUM(1]  OlOOO  THEN  STATSONATTS; 

FOR  I ;=  1  TO  10  DO  VJRITELN; 

WRITELN ( •  INDIVIDUAL  SCORES ’ ) ; 

WRITELNC  - •); 

FOR  I;«l  TO  3  DO  WRITELN; 

WRITE ( 'LNR' ) ; 

FOR  I:*l  TO  NUMOFTALKERS  DO  WRITEC  T’,I:1); 

WRITELN; 

FOR  I;»l  TO  NUMOFLISNERS  DO 
BEGIN 

WRITE(I;2,’  '); 

FOR  J;=l  TO  NUMOFTALKERS  DO  WRITE(SINREPMAT[ I, J) : 6; 1) 
WRITELN; 

END; 

WRITELN; 

WRITELN ('LNR  MEAN  OVER  TKR ; 

FOR  I:-l  TO  NUMOFLISNERS  DO 

WRITELN(I}2,LMEANS(I):13:1); 

WRITELN; 

FOR  ls«l  TO  S  DO  WRITELN; 

WRITEI 'TKR' ) ; 

FOR  I:«l  TO  NUMOFLISNERS  DO 
WRITEC  L',I:1); 

WRITELN; 

FOR  I:»l  TO  NUMOFTALKERS  DO 
BEGIN 

WRITE (I: 2,'  '); 

FOR  J;«l  TO  NUMOFLISNERS  DO  WRITE{SINREPMAT{ J, I ) ; 6: 1) 
WRITELN; 

END; 

WRITELN; 

WRITELN ('TKR  MEAN  OVER  LNR'); 

FOR  I;-l  TO  NUMOFTALKERS  DO 

WRITELN(I:2,TMEANSII) :13;1); 

WRITELN; 

GRANDMEAN:»0; 

FOR  Jj«l  TO  NUMOFTALKERS  DO 

GRANDMEAN:»GRANDHEAN  +  TMEANS(J); 

GRANDHEAN ; -GRANDHEAN/NUMOFTALKERS ; 

FOR  Ij-1  TO  5  DO  WRITELN; 

WRITELN (’DRT  SCORE  «  ' ,GRANDMEAN: 5; 1) ; 

FOR  I:»l  TO  5  DO  WRITELN; 

END; 


(*  Main  Program  *) 

BEGIN 

WRITELN; 

FOR  I:=l  TO  3  DO  WRITELN; 
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INITIALISE; 

READSEQFILE; 

READRESFIL; 

SINGLEREPS; 

CALCLMEANS; 

CALCTMEANS; 

IF  (NUMOFTALKERS  >  2)  AND 
(NUMOFLISNERS  >  2) 

THEN  ANOVAR 

ELSE  BASIC_AND_ATT_STATS; 
UPDATECONDMEANS ; 

END. 
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PROGRAM  DRTFILE(INPUT, OUTPUT, AFL,BFL,DRTTMPFLl,DRTEDA,CONDMEANS) ; 


(*  Input:  Data  from  the  file  Condmeans  created  *) 

(*  by  DRTSTAT.  *) 

(*  Processing:  Extract  required  results,  transpose  *) 

(*  into  chars,  reformat  and  store  for  *) 

(*  subsequent  use  by  PANOVUN.  *) 

CONST 

CREWNAMES»100; 

DEPTH0FNESTING>=15  ; 

MAXDATASTRG«=100; 

MAXTITLESTRG=100 ; 


TYPE 

ARROFMEANS*ARRAY ( 1 . . 12 )  OF  REAL; 

STROFTWO=PACKED  ARRAY (1.. 2)  OF  CHAR; 

CONDTYP-PACKED  ARRAY (1.. 4)  OF  CHAR; 

ATTYP=PACKED  ARRAY (1.. 3)  OF  CHAR; 

REALCHRS=PACKED  ARRAY (1.. 5]  OF  CHAR; 

NAHTYP*PACKED  ARRAY {1.. 6)  OF  CHAR; 

STROFTHREE-PACKED  ARRAY (1.. 3]  OF  CHAR; 

CONDREC=RECORD 

CONDITION; (KEY (0)1  PACKED  ARRAYtl..71  OF  CHAR; 
LNAMS : ARRAY (1.. 12)  OF  CONDTYP; 

LISNEROAMS : ARROFMEANS ; 

ATTRIBUTES:ARRAY(1..12,1..2,1..6)  OF  REAL; 
TLKRESMAT:ARRAY(1..12,1..12)  OF  REAL; 

TNAMS ; ARRAY { 1 .. 12 )  OF  STROFTHREE; 

TLKNUM: INTEGER; 

END; 


VAR 

CONDHEANS;FILE  OF  CONDREC; 

AFL , BFL , DRTTHPFLl , DRTEDA : TEXT ; 

CONDLEN, ATTLEN , LNUM, I , NUH: INTEGER; 
CONDNUMSTRG,ATTSTRG; VARYING (1000)  OF  CHAR; 

COND; CONDREC; 

PAM:ARRAY(1..12,1..3,1..6)  OF  REAL; 

CNUM; CONDTYP; 

CH;CHAR; 

FLSNOTFOUND;INTEGER; 

NAMES, PERSONS: ARRAY (1..CREWNAMES)  OF  NAMTYP; 

CREW, COL, EDACOLS; INTEGER; 

NESTING, NESTKAX: ARRAY (1..DEPTHOFNESTING)  OF  INTEGER; 
LISNERMISSING: ARRAY (1..CREHNAHES)  OF  INTEGER; 

INDEX, LIDNUH: ARRAY (1..CREWNAMES)  OF  STROFTWO; 
REQUDLNRS : ARRAY ( 1. .CREWNAMES)  OF  NAMTYP; 
EXPECTED_VALS__PER_LNR;  INTEGER; 

PERCENT; REAL; 

FUNCTION  INTTOCHR(INTIN:INTEGER) ;STR0FTW0; 

{*  turns  integers {<99)  into  chars  *) 

VAR 
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ICHRStARRAYll. .2}  OF  CHAR; 

ILOCSTRG : STROFTWO ; 

BEGIN 

ICHRS[1]:«CHR(48+(INTIN  DIV  10)); 

ICHRS(2):*CHR(48+(INT1N  REM  10)); 

IF  ICHRS(1]«’0’  THEN  ICHRS(1]:**  •; 

PACK ( ICHRS , 1 , ILOCSTRG) ; 

INTTOCHR : ^ILOCSTRG ; 

END; 

FUNCTION  CHRTOINT( CHIN: STROFTWO) JINTEGER; 

(*  turns  2  char  string*  into  integers  *) 

VAR 

INTOUT: INTEGER; 

CHRINTS:ARRAy(l. .2]  OF  CHAR; 

BEGIN 

UNPACK(CHIN,CHRINTS, 1) ; 

IF  CHRINTS(2]='  ' 

THEN 

BEGIN  CHRINTS(2J:=CHRINTS(1);CHRINTS(1);*='0'  END; 
INTOUT:-{10*(ORD(CHRINTS(l))-48))  +  (ORD(CHRINTSI2))-48); 
CHRTOINT:«INTOUT 

END; 

FUNCTION  TOCHARS(LSNR;INTEGER;R:REAL) :REALCHRS; 

(*  turns  poB  or  neg  real  numbers  i  to  chars  of  the  form  SAA.A  *) 
VAR 

ICER: INTEGER; 

CHRS:ARRAY(1. .5)  OF  CHAR; 

LOCALSTR:REALCHRS; 

BEGIN 

R:-(ROUND(R*10) )/10; 

IF  R-lOO  THEN  TOCHARS:-'  100  * 

ELSE 

IF  (R-O)  AND  (COND.LNAMS(LSNRJ-' - ’) 

THEN  TOCHARS:-’  -999* 

ELSE 

BEGIN 

IF  R<0  THEN 

BEGIN 

CHRSilJi*'-’; 

R:-R*{-1) 

END 

ELSE  CHRS(1J:«*  ’; 

IGER:-TRUNC(R*10); 

CHRS(2]:-CHR(48+(ICER  DIV  100)); 
CHRS[3):-CHR((48+(ICER  DIV  10)  REM  10)); 
CHRSt4):-'. '; 

CHRS(5):»CHR(48+(IGER  REM  10)>; 

IF  CHRS[2J»'0'  THEN  CHRS(2):»*  '; 

IF  (CHRS(1J='-*)  AND  {CHRS[2]='  ') 

THEN 

BEGIN 

CHRSll):®'  •; 
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CHRS[2] 

END; 

PACK { CHRS , 1 , LOCALSTR ) ; 

TOCHARS ; =LOCALSTR 
END; 

END; 

PRCXJEDURE  ENTERNAMES; 

(*  read  listener  names  at  start  of  each  condition  *) 

VAR 

I , TIMESTHRU ; INTEGER; 

MATCHED : BOOLEAN ; 

BEGIN 

FOR  I;»l  12  DO 
BEGIN 

TIMESTHRU  :•=!; 

READ ( DRTTMPFLl , NAMES ( I ) ) ; 

MATCHED: -FALSE; 

IF  NAMES[I1<>'  - • 

THEN 

BEGIN 

WHILE  (NOT  MATCHED)  AND  (TIMESTHRU<-CREW)  DO 
IF  NAMES! I]  «  PERSONS (TIMESTHRU) 

THEN 

BEGIN 

MATCHED: -TRUE; 

INDEX ( 1 1 : -LIDNUH ( TIMESTHRU ) 

END 

ELSE 

TIMESTHRU : -TIMESTHRU+1 ; 

IF  NOT  MATCHED 

THEN  WRITELN( 'ERROR,  Unrecognised  listener 

•,NAMES(I)); 

END; 

END; 

READLH (DRTTMPFLl ) ; 

END; 


PROCEDURE  UPDATE(VAR  CURRENT, MAX  :  INTEGER); 
BEGIN 

IF  CURRENT<MAX 
THEN 

CURRENT: -CURRENT+1 

ELSE 

BEGIN 

CURRENT:-!; 

COL:-COL+l; 

UPDATE (NESTING (COL) ,NESTMAX(COL) ) ; 
END; 

END; 


name 
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PRCXJEDURE  OUTPTNESTING; 


BEGIN 

COL;*l; 

UPDATE ( NESTING { COL  J , NESTMAX [ COL ] ) ; 

FOR  I:-l  TO  EDACOLS+1  DO  TOITE{BFL, ’  * ,INTTOCHR(NESTINGtI) ) ) ; 

END; 


PROCEDURE  TOEDAFORM; 

(*  rsfonnat  temp  file  to  EDA  format  *) 

VAR 

X,Y,I,J,K,L, CONDSETS , NUM , REQLNR , TEMPI : INTEGER ; 

DATANAMES: VARYING [MAXDATASTRG]  OF  CHAR; 

TITLE: VARYING (MAXTITLESTRG)  OF  CHAR; 

DRTiNAMTYP; 

SCORE;ARRAY(l. .12)  OF  NAMTYP; 

WHOBY : ARRAY (1. .12)  OF  INTEGER; 

CH:CHAR; 

OK:BOOLEAN; 

REQULNRNUMS;ARRAY(1. .12)  OF  STROFTWO; 

TEMP2; NAMTYP; 

BEGIN 

(*  read  control  file  •) 

I:«l  TO  10  DO  READ{AFL,CH); 

P-SADLN  { AFL ,  REQLNR ) ; 

FOR  I:«l  TO  REQLNR  DO 
BEGIN 

READLN ( AFL , REQUDLNRS ( I ) ) ; 

LISNERMISSING(I):«0; 

WHOBY(I) !-0 
END; 

READLN(AFL,CH); 

IF  CH  <> 

THEN  WRITELN  ('ERROR  **  Caused  by  incorrect  number  of  listener 

IDS'); 

READLN ( AFL); 

READLN ( AFL , EDACOLS ) ; 

READLN (AFL, DATANAMES ) ; 

NESTMAX(l)  :»  REQLNR; 

FOR  I:»2  TO  EDACOLS+1  DO 
READ (AFL , NESTMAX ( I ) ) ; 

READLN (AFL); 

READLN (AFL, TITLE); 

READLN (AFL, CH); 

IF  CHO'*'  THEN  vmiTELN( 'ERROR  **  Caused  by  too  many  format 
lines' ) ; 

REP7JAT 

READLN (AFL, CH) 

UNTIL  CH*'*'; 

FOR  I:=l  TO  16  DO  READ(AFL,CH) ; 

READLN (AFL, CREW) ; 
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FOR  I:*l  TO  CREW  DO 

READLN(AFL,PERS0NS(I],CH,L1DNUM[I]);  (*  read  listener  ids  *) 
FOR  I:=l  TO  KEQLNR  DO  (*  note  requd  listener 

ids  *) 

FOR  J;*l  TO  CREW  DO  IF  REQUDLNRS ( I ] -PERSONS [ J ]  THEN 

REQULNRNUMS [ I ] : =LIDNUM ( J ] ; 

NESTING! 1] ;-0; 

FOR  I:«2  TO  EDACOLS+1  DO  NESTING ( I ]: «1 ; 

WRITELN(BFL,INTTOCHR{EDACOLS+3));  {*  write  to  opfile  *) 

WRITELN(BFL, ' :"DRT  ","INDIV  ", "POSIT  " , ' , DATANAMES ) ; 

WRITELN ( BFL , TITLE ) ; 

RESET (DRTTMPFLl) ; 

CONDSETS;-!; 

FOR  I: -3  TO  EDACOLS+1  DO 

CONDSETS ; -CONDSETS"NESTMAX ( I 1 ; 

EXPECTED__VALS_PER_LNR: -NESTMAX ( 2 ) *CONDSETS ; 

FOR  It-l  TO  CONDSETS  DO 
BEGIN 

ENTERNAMES ; 

FOR  J;-l  TO  NESTMAX(2)  DO 
BEGIN 
NUM?=1; 

FOR  K:=l  TO  12  DO  WHOBy(K)!=0; 

FOR  K;-l  TO  12  DO 
BEGIN 

READ ( DRTTMPFLl , DRT ) ; 

OK : -FALSE; 

FOR  Ls-1  TO  REQLNR  DO  IF  NAMES ( K ) -REQUDLNRS ( L J 

THEN  OK : -TRUE; 

IF  OK 
THEN 

BEGIN 

(*  build  up  req  listeners  for  sorting  *) 
SCORE[NUHl:-DRT; 

WHOBY(NUM):»CHRTOINT(INDEX(K) ) ; 

NUM:*NUH+1 

END; 


(*  check  each  req  listener  has  a  score  *) 

WHILE  (NUM-1)<REQLNR  DO 
FOR  X:«l  TO  REQLNR  DO 
BEGIN 
OK ; -FALSE; 

FOR  yj-1  TO  (NUM-1)  DO 

IF  WHOBY{YJ«CHRTOINT{ REQULNRNUMS (X)) 
THEN  OK : -TRUE; 

IF  NOT  OK 
THEN 

BEGIN 

WHOBY(NUMl ;»CHRTOINT (REQULNRNUMS (X] ) ; 
SCORE(NUM):»'  -999’; 

LISNERMISSING (X ) : -LISNERMISSING ( X ) +1 ; 

NUM:=HUH+1 

END; 
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END; 

(*  sort  *) 

FOR  K;=l  TO  REQLNR-1  DO 

FOR  L;=l  TO  REQLNR-1  DO 

IF  WHOBY(L]>WHOBYIL+l} 

THEN 

BEGIN 

TEMPl;*WHOBY(L); 

WHOBir  [  L  ] :  «WHOBY  I L+1  ]  ; 

WHOBY ( L+1 ] : -TEMPI ; 

TEMP2 ; -SCORE (L); 

SCORE [L] ;-SCORE(L+l); 

SCORE ( L+1 ]:-TEMP2 
END; 

(*  print  *) 

FOR  K!=1  to  REQLNR  DO 
BEGIN 

V?RITE(BFL,SC0RE(K),  •  •  ,INTTOCHR(WHOBY(K) ) )  ; 

OUTPTNESTING; 

WRITELN(BFL) 

END; 

READLN ( DRTTMPFLl ) ; 

END; 

READLN ( DRTTMPFLl ) ; 

END; 

FOR  I:»l  TO  EDACOLS+3  DO  WRITE ( BFL, '-1  '); 

WRITELN(BFL) ; 

WRITELN(BFL, '-1' ); 

END; 

PROCEDURE  DOEDAINPT; 

(*  create  input  file  for  panovun  *) 

VAR 

I : INTEGER; 

CHrCHAR; 

BEGIN 

REWRITE (DRTEDA); 

RESET (AFL); 

FOR  It-1  TO  5  DO 
BEGIN 
CH;»’ 

WHILE  CHO’*'  DO  READ(AFL,CH); 

READLN (AFL) 

END; 

READLN (AFL); 

CH;»'  •; 

WHILE  CHO'*'  DO 
BEGIN 

WHILE  NOT  EOLN(AFL)  DO 
BEGIN 

READ(AFL,CH); 

IF  CHO'*’  THEN  WRITE  (DRTEDA,  CH) 
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END; 

READLN{AFL) ; 
WRITELN(DRTEDA) ; 
END; 

END; 


PROCEDURE  OUTDATA(ATT:ATTyP); 

(*  decode  attribute  required  and  output  *) 


VAR 


L0CALARR:ARRAY[1. .3]  OF  CHAR; 

TLKPOS , TPOS , TABNO, I , J : INTEGER; 

PAORM:CHAR; 

BEGIN 

UNPACK ( ATT , LOCALARR , 1 ) ; 

PAORM ; =LOCALARR (3); 

TABNO:*ORD(LOCALARR(l) )+ORD(LOCALARR[2) ) ; 
IF  PAORM=’P’  THEN  J:=l 

ELSE  IF  PAORM=’A’  THEN  J:«2 
ELSE  J;=3; 


WITH  COND  DO 

CASE  TABNO  OF 

144:  BEGIN 

FOR  I:=l  TO  12  DO 

WRITE (DRTTMPFLl , ’  ’ , TOCHARS ( I , LISNEROAMS ( I ) ) ) ; 

END; 

165:  BEGIN 

FOR  I:«=l  TO  12  DO 

WRITE (DRTTMPFLl, ’  • ,TOCHARS(I,PAM(I, J, 1) ) ) ; 

END; 


143: 


168: 


156: 


153: 


146: 


159: 


BEGIN 

FOR  I:»l  TO  12  DO 

WRITE (DRTTMPFLl, ’  • ,TOCHARS ( I , PAMJ I , J, 2 ) ) > ; 

END; 


BEGIN 

FOR  I:=l  TO  12  DO 

WRITE ( DRTTMPFLl , ’  ' , TOCHARS ( I , PAM ( I , J , 3 ] ) ) ; 

END; 

BEGIN 

FOR  I:=l  TO  12  DO 

WRITE (DRTTMPFLl, '  * ,TOCHARS(I,PAM(I, J,4) ) ) ; 

END; 

BEGIN 

FOR  I:«l  TO  12  DO 

WRITE (DRTTMPFLl, '  * ,TOCHARS(I,PAM(I,J,5) ) ) ; 

END; 

BEGIN 

FOR  I;*l  TO  12  DO 

WRITE (DRTTMPFLl, *  * ,TOCHARS ( I , PAMt I , J, 6 ) ) ) ; 

END; 

BEGIN 

FOR  I:=l  TO  TLKNUM  DO 
BEGIN 
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FOR  J:=l  TO  12  DO 

• ,TOCHARS(J,TLKRESMAT{J,I) ) ) ; 

WRITELN(DRTTMPFLl) ; 
END; 

END; 


WRITE  (DRTTMPFLl,  ’ 


END; 

IF  TABNO  <>  159  THEN  WRITELN ( DRTTMPFLl ) ; 

END; 

PROCEDURE  SELECTDATA  ( COND  ;'C0NDREC )  ; 

(*  write  heading  for  block  of  data  and  loop  thru  attributes  *) 
(*  calling  OUTDATA.  *) 

VAR 

ANUM: INTEGER; 

ATTREQ;ATTyP; 

BEGIN 

FOR  I:=l  TO  12  DO 

WRITE (DRTTMPFLl, ’  ' , COND . LNAMS ( I ] ); 

WRITELN ( DRTTMPFLl ) ; 

ANUM: =4; 

REPEAT 

ATTREQ :  >=SUBSTR  ( ATTSTRG ,  ANUM-3 ,3); 

OUTDATA (ATTREQ) ; 

ANUM:=ANUM+4 
UNTIL  ANUM  >  ATTLEN; 

WRITELN (DRTTMPFLl) 

END; 

PROCEDURE  RETRI EVEREC ( CONDNUMB : CONDTY P ) ; 

(*  locate  required  conditions  using  EXPID  as  key,  generate  *} 
(*  means  from  Present  and  Absent,  and  call  SELECTDATA  •) 

VAR 

EXPID: PACKED  ARRAY (1., 7)  OF  CHAR; 

I,  J,  K  .‘INTEGER; 

BEGIN 

EXPID : * ' CTL ' +CONDNUHB ; 

COND . CONDITION : »EXPID ; 

OPEN(CONDMEANS, 

HISTORY :=UNKNOWN, 

ORGANIZATION : -INDEXED , 

ACCESS_METHOD; -KEYED) ; 

RESETK ( CONDMEANS , 0 ) ; 

FINDK ( CONDMEANS , 0 , COND . CONDITION ) ; 

IF  NOT  UFB( CONDMEANS) 

THEN 

BEGIN 

READ ( CONDMEANS , COND ) ; 

FOR  I:-l  TO  12  DO 

FOR  K:»l  TO  6  DO 
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BEGIN 

FOR  J:=l  TO  2  DO 

PAM[I,J,K]:=COND.ATTRIBUTES[I,J,K]; 

PAM[I,3,K]:=(PAM(I,l,K)+PAM[I,2,K])/2; 

END; 

SELECTDATA ( COND ) ; 

END 

ELSE 

BEGIN 

WRITSLN(EXPID, •**  not  found  **'); 

PLSNOTFOUND ; -rLSNOTFOUND+l ; 

WRITELN (DRTTMPFLl ) 

END; 

CLOSE ( CONDMEANS ) ; 

END; 


PROCEDURE  DECODESTRGS; 

(*  decode  string  of  conditions  requd,  and  loop  thru  calling  RETRIEVEREC 

*) 

VAR 

NUMl , NUM2 , I : INTEGER; 

DELIM! PACKED  ARRAY (1..1]  OF  CHAR; 

CONDARR,LASTARR:ARRAY(l. .4)  OF  CHAR; 

ATTARR:ARRAy(l. .2]  OF  CHAR; 

CONVS , CONDNUM , LASTNUM : CONDTYP ; 

BEGIN 

CONDLEN : -LENGTH { CONDNUMSTRG ) ; 

ATTLEN : -LENGTH ( ATTSTRG ) ; 

LNUH:=1; 

REPEAT 

CONDNUM : -SUBSTR ( CONDNUMSTRG , LNUM , 4 ) ; 

DELIM : -SUBSTR ( CONDNUMSTRG , LNUM+4 , 1 ) ; 

IF  DELIM-'-' 

THEN 

BEGIN 

LASTNUM : -SUBSTR ( CONDNUMSTRG , LNUM+5 , 4 ) ; 
UNFACK(CONDNUH,CONDARR, 1) ; 

UNPACK ( LASTNUM , LASTARR, 1 ) ; 

NUMl : -1000+100* (ORD (CONDARR [ 2 ) ) -48 ) +10* (ORD ( CONDARR ( 3 ) ) -48 ) 

+ORD ( CONDARR [ 4 ) ) -48 ; 

NUM2 t -1000+100* (ORD (LASTARR( 2 1 ) -48) +10* (ORD (LASTARR! 3 )) -48 ) 

+ORD ( LASTARR ( 4 ] ) -4 8 ; 

FOR  I: -NUMl  TO  H0H2  DO 
BEGIN 

CONDARR[4)z»CHR(48+(I  REM  10)); 
CONDARR(3];«CHR(48+((I  REM  100)  DIV  10)); 
CONDARR(2);«CHR(48+((I  REM  1000)  DIV  100)); 

PACK ( CONDARR , 1 , CONDNUM ) ; 

RETRIEVEREC (CONDNUM) ; 
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END? 

LNUM:=LNUM+5; 

END 

ELSE 

RETRIEVEREC (CONDNUM) ; 

I,NUM;»LNUM+5 

UNTIL  LNUM  >  COKDLEN; 

END; 

BEGIN 

OPEN(AFL, • infile', 

OLD)  ; 

RESET (AFL); 

OPEN(BFL, 'outfile'); 

REWRITE ( BFL); 

REWRITE (DRTTMPFLl ) ; 

CONDNUMSTRG !  =  "  ; ATTSTRG :  =  ' ' ; FLSNOTFOUND : =0 ; 

FOR  I:=l  TO  3  DO  READLN(AFL); 

READ (AFL, CH) ; 

WHILE  CHO'*'  DO 
BEGIN 

IF  ORD{CH)=13  THEN  CH:=','; 

ATTSTRG ; »ATTSTRG+CH ; 

READ (AFL, CH) ; 

END; 

READLN { AFL ) ; READLN ( AFL ) ; 

READ(AFL,CH) ; 

WHILE  CH<>'»'  DO 
BEGIN 

IF  ORD(CH)«13  THEN  CH:=','; 

CONDNUMSTRG : =CONDNUMSTRG+CH ; 

READ(AFL,CH) ; 

END; 

READLN (AFL); 

FOR  Ij«=l  TO  4  DO  WRITELN; 

WRITELNCOK,  running  DRTFILE'); 

WRITELN; 

DECODESTRGS; 

IF  FLSNOTFOUND>0 
THEN 

BEGIN 

WRITELN; 

WRITELN ( 'Error J  **', FLSNOTFOUND ; 3 , '  files  not  found  **'); 
WRITELN; 

END 

ELSE 

BEGIN 

WRITELNCOK,  all  records  found,  creating  output  file'); 
WRITELN; 

TOEDAFORM; 

WRITELNCOK,  output  file  created'); 

WRITELN; 

FOR  I:>=1  TO  CREW  DO 

IF  LISNERMISSING(I)>0 
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THEN 

BEGIN 

PERCENT;=(L1SNERMISSING[I]  /  EXPECTED_VALS_PER_LNR) *100; 

WRITELN{ 'Warning  , * , REQUDLNRS ( I ] , •  has 

LISNERMISSING{I]53,'  (',PERCENT;3:1,'%)  missing 

values' ) ; 

END; 

DOEDAINPT; 

END; 


PROGRAM  DRTEDIT ( INPUT, OUTPUT, CONDKSANS ) ; 


TYPE 

CONDTYP=PACKED  ARRAY [1.. 4]  OF  CHAR; 

CONDNAMETYP=PACKED  ARRAY (1.. 7 J  OF  CHAR; 
S,TROFTHREE»=PACKSD  ARRAY [1.. 3)  OF  CHAR; 
ARROFMEANS*ARRAY(l. il2]  OF  REAL; 

CONDREC«RECORD 

CONDITION: (KEY (0) ]  PACKED  ARRAYfl..?]  OF  CHAR; 
LNAMS:ARRAY(1. .12]  OF  CONDTYP; 

LISNEP.OAMS :  ARROFKEANS ; 

ATTRIBUTES;ARRAY[1..12,1..2,1..6}  OF  REAL; 
TLKRESMAT : ARRAY (1.. 12,1.. 12]  OF  REAL; 

TNAMS : ARRAY (1.. 12]  OF  STROFTHREE; 

TLKNUM: INTEGER; 

END; 


VAR 

CONDMEANS:FILE  OF  CONDREC; 

CHOICE: INTEGER; 

COND: CONDREC; 

PROCEDURE  DISPINDEX; 

VAR 

WRI TS; INTEGER ; 

NUM: CONDTYP; 

BEGIN 

WRITELN; 

RESET (CONDMEANS); 

WRITS:*-1; 

WHILE  NOT  EOF (CONDMEANS)  DO 
BEGIN 

READ ( CONDMEANS , COND ) ; 

NUM : -SUBSTR ( COND . CONDITION ,4,4); 
WRITS :»WRITS+1; 

IF  WRITS>12 
THEN 

BEGIN 

WRITELN; 

WRITS :»0 
END; 

WRITE (NUM, •  '); 

END; 

WRITELN; 

CLOSE (CONDMEANS); 

END; 

PROCEDURE  DELENTRY; 

VAR 

CNUM: CONDTYP; 

EXPID : CONDNAMETYP ; 

BEGIN 

OPEN (CONDMEANS, 

HISTORY : =UNKNOWN , 
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ORGANIZATION :=INDEXED, 

ACCESS_METHOD ; *KEYED ) ; 

WRITE ('Enter  Condition  Number  to  be  Deleted:  '); 
READLN(CNUM); 

EXPID;*'CTL'+  CNUM; 

FINDK ( CONDMEANS , 0 , EXPID ) ; 

IF  NOT  UFB( CONDMEANS) 

THEN 

BEGIN 

DELETE ( CONDMEANS ) ; 

WRITELN (EXPID, •  deleted  OK') 

END 

ELSE 

WRITELN ( 'Condition  ', EXPID, '  not  on  file'); 
CLOSE (CONDMEANS) 

END; 

PROCEDURE  DISPMENU; 

BEGIN 

WRITELN; 

WRITELN ( 'DRTEDIT' ) ; 

WRITELN  ( ' - ' ) ; 

WRITELN; 

WRITELN ('1.  Display  Index'); 

WRITELN('2.  Delete  an  entry'); 

WRITELN ( ' 3 .  End ' ) ; 

WRITELN; 

WRITE ( 'Select. .  '); 

END; 

BEGIN 

CHOICE :»0; 

WHILE  CH0ICEO3  DO 
BEGIN 
DISPMENU; 

READLN ( CHOICE ) ; 

CASE  CHOICE  OF 
1:  DISPINDEX; 

2:  DELENTRY; 

END; 

END; 

END. 
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A3  Sample  Program  Outputs 
Output  from  DRTI.COM 
DRT  STATISTICS  PROGRAM 


Enter  filename;  Enter  filename; 


Control  Pile;  CTL3139 


Source:  Multi-Listener  Master  V4 
Date:  9 -NOV-8 9 

Input  Device;  BSK  WB  TAPE 
Input  Noise:  NONE 

Link  Device:  SPECIAL  CODER 

Output  Device;  BT  HANDSET 
Output  Noiqe:  6SDB(A)USASI 

•  Optional  Header  Text 
HANDSET  TEST 
TAPE  8 

Sequence  Type:  DYNA 
Filler  Display;  Enabled 

Words  File:  DYNAWRD 

Sequences  File;  DYNASEQ 

DRTs:  10 


I  Talker 

List 

RLP 

306A.1 

RLP 

306A.2 

RB 

308A.1 

RB 

308A.2 

KCL 

315A.1 

MCL 

31SA.2 

PRW 

312a. 1 

PRW 

312A.2 

GAP 

314A.1 

GAP 

314A.2 

Listeners: 

10 

1  Listener 

Station 

ABBOTT 

1 

BERRY 

2 

BLACK 

3 

MCINTOSH 

4 

RUSBY 

S 

HOOK 

6 

JONl 

7 

TAYLOR 

8 

JON2 

9 

SLATER 

10 

CND3139 
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TUKEY  TEST  FOR  NONADDATIVITY 


Source 

Sura  of 

Degrees  of 

Mean 

Squares 

Freedom 

Square 

LISNERS 

477.09 

9 

53.01 

TALKERS 

624.88 

4 

156.22 

RESIDUAL 

460.19 

36 

NONADD 

3.38 

1 

3.38 

BALANCE 

456.81 

35 

13.05 

ANALYSIS  OF  VARIANCE 


Source 

Sum  of 

Degrees  of 

Mean 

F-Ratio 

Squares 

Freedom 

Square 

LISNERS 

954.25 

9 

106.03 

1.99 

TALKERS 

1249.75 

4 

312.44 

5.87 

LN  •  TK 

920.38 

36 

25.57 

0.48 

ERROR 

2660.56 

50 

53.21 

TOTAL 

5784.94 

99 

NEWMAN  -  KEULS  TEST 


LISTENERS  : 

238461  10  975 


TALKERS  : 

4  3  5  1  2 


ATTRIBUTE  SCORES 


Present  Absent  Mean 


mean 

SE 

mean 

SE 

mean 

SE 

Voicing 

87.5 

2.5 

70.5 

4.7 

79.0 

3.2 

Nasality 

96.5 

0.6 

84.3 

3.5 

90.4 

2.2 

Sustention 

72.8 

2.8 

69.3 

2.8 

71.0 

2.0 

Sibilation 

82.5 

1.9 

81.5 

2.2 

82.0 

1.4 

Graveness 

61.0 

2.4 

62.3 

2.8 

61.6 

1.8 

Compactness 

90.0 

2.6 

80.5 

3.9 

85.3 

2.5 

INDIVIDUAL  SCORES 


LNR 

T1 

T2 

T3 

T4 

T5 

SD 

1 

79.2 

82.3 

72.9 

76.0 

86.5 

5.3 

2 

76.0 

80.2 

67.7 

61.5 

71.9 

7.3 

3 

82.3 

78.1 

72.9 

70.8 

71.9 

4.9 

4 

77.1 

81.3 

76.0 

80.2 

76.0 

2.4 

5 

84.4 

86.5 

76.0 

78.1 

86.5 

4.9 

6 

83.3 

82.3 

76.0 

77.1 

78.1 

3.2 

7 

84.4 

85.4 

77.1 

68.8 

88.5 

8.0 

8 

75.0 

80.2 

79.2 

65.6 

77.1 

5.8 

9 

83.3 

81.3 

76.0 

76.0 

85.4 

4.3 

10 

81.3 

83.3 

76.0 

76.0 

81.3 

3.3 

LNR  MEAN  OVER  TKR 

1  79.4 

2  71.5 

3  75.2 

4  78.1 

5  82.3 

6  79.4 

7  80.8 

8  75.4 

9  80.4 

10  79.6 

STD.  DEVN.  =  3.3 

STD.  ERROR  =  1.0 
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TKR  LI  L2  L3  L4  L5 

1  79.2  76.0  82.3  77.1  84.4 

2  82.3  80.2  78.1  81.3  86.5 

3  72.9  67.7  72.9  76.0  76.0 

4  76.0  61.5  70.8  80.2  78.1 

5  86.5  71.9  71.9  76.0  86.5 

TKR  MEAN  OVER  LNR 
1  80.6 
2  82.1 

3  75.0 

4  73.0 

5  80.3 

STD.  DEVN.  >=  4.0 

STD.  ERROR  =  1.8 


DRT  SCORE  =78.2 
STANDARD  ERROR  =  1.0 


L6  L7  L8  L9  LIO  SD 

83.3  84.4  75.0  83.3  81.3  3.5 

82.3  85.4  80.2  81.3  83.3  2.5 

76.0  77.1  79.2  76.0  76.0  3.1 

77.1  68.8  65.6  76.0  76.0  6.1 

78.1  88.5  77.1  85.4  81.3  6.2 
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Output  from  DRT2.COM  (written  to  DRTTMPFL2.DAT) 
5 

:"DRT  ","INDIV  ", "POSIT  ", "ATTR", "COND" 

TITLE 

62.5  14  1  1  1 

75.0  14  1  2  1 

75.0  14  1  3  1 

etc. . . 

75.0  15  2  4  4 

46.9  15  2  5  4 

75.0  15  2  6  4 

~1  -1  -1  -1  -1 

-1 
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Output  from  DRT3.COM 


CCF>  edrt3 
DRTEDIT 


1.  Display  Index 

2.  Delete  an  entry 

3.  End 


Select..  1 

1191 

1192 

1193 

1194 

1195 

1196 

1197 

1199 

1200 

1201 

1202 

1203 

1204 

1205 

1206 

1207 

1208 

1219 

1220 

1221 

1222 

1223 

1224 

1225 

1226 

1227 

1228 

1229 

1230 

1231 

1232 

1233 

1234 

1235 

1236 

1237 

1238 

1239 

1240 

1241 

1242 

1243 

1244 

1245 

1246 

1247 

1248 

1249 

1250 

1251 

1252 

1253 

1254 

1255 

1256 

1257 

1258 

1259 

1260 

1261 

1262 

1263 

1264 

1265 

1266 

DRTEDIT 


1.  Display  Index 

2.  Delete  an  entry 

3.  End 


Select..  3 


A43 


A4  Sample  Control  File  for  DRT2 


Example  Control  File  (CF....) 

1  This  file  is  for  EXAMPLE  purposes  only,  and  will  not  run. 

1  Lines  starting  with  1  are  used  for  comment  here,  but  will  cause  fatal 
errors  at  runtime. 

1 

1  Make  the  control  file  n2une  the  same  as  the  filename  used  to  store 
this  file 
1 

1  N.B.lil  The  “star"  symbols  must  seperate  the  different  entries. 

! 

Control  File;  CFwhatever 
* 

1 

1  These  numbers  after  the  colon  must  tally  with  the  number  of  lines 
below  them. 

1 

Attributes;! 

TKR 

* 

J  List  CND  numbers  here.  These  are  in  the  order  corresponding  to  the 
1  Design  statement  used  in  PANOVUN. 

1 

Conditions: 24 

1338 

1349 

1342 

1353 

1335 
1348 
1334 

1345 

1331 
1344 

1339 
1352 
1337 

1350 
1341 

1354 

1336 
1347 
1333 

1346 

1332 

1343 

1340 

1351 
* 

I 

t  Select  listeners  whose  results  you  wish  to  use. 

I  Some  degree  of  "  missing  values'*  (less  than  10%)  is  acceptable 
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1  you  MUST  LEAVE  TWO  SPACES  in  front  of  the  four  letter  name. 

1 

Listeners: 12 
CHAM 
PRAT 
WHEA 
SLAT 
PAWS 
BERK 
FERR 
PERK 
SIMO 
KNUT 
TAYL 
PLOW 

* 

1 

1  First  line  is  ALWAYS  "Format; 4" 

1  Second  line  is  the  number  of  variables, 
i  Third  line  is  the  nesting  of  the  variables. 

I  Fourth  line  is  the  quantity  of  each  of  the  variables. 

1  Fith  line  is  a  title  that  will  appear  on  the  printout. 

1 

Format ; 4 
5 

"TALKER", "REPLIC", "UTTER  ", "CODER  ","SNR 
5  2  2  3  2 
Title  line 

A 

I 

!  These  statements  drive  PANOVUN,  an  analysis  of  variance  program 
1  Written  by  the  Institute  of  Aviation  medecine. 

I 

Anovar  Design: - 
Another  title 
Y 
1 
H 
:E 
1 
E 

D  Put  your  design  statement  here 
A(IY,V— 999)DRT 
E 

Q 

* 

I 

1  Listener  ident.  numbers  are  put  here. 

1  Last  listener  is  always  " - " 

1  Further  listeners  may  be  added  as  required. 

I  but  leave  TWO  SPACES 
1 
1 
I 


A45 


Listener  Idents:14 
CHAM  1 
PRAT  2 
VfHEA  3 
SLAT  4 
PAWS  5 
BERK  6 
KERR  7 
PERK  8 
SIMO  9 
KNUT  10 
TAYL  11 
PLOW  12 
ROBE  13 
- 14 
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